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REGRESSION ANALYSIS PROCEDURES FOR THE EVALUATION 
OF TRACKING SYSTEM MEASUREMENT ERRORS 

SUMMARY 


The TEMS Multiple Regression Analysis Method for evaluating systematic 
errors in measurements obtained from various tracking systems is presented. 
The mathematical procedures in the method include a rigorous least squares 
adjustment of error model parameters with constraints. The mathematics for 
using a priori values for these parameters and their variances is also in- 
cluded. Truncated tracker error models for representing the systematic 
errors are established using the TEMS method in conjunction with a stepwise 
regression procedure. The basic approach in the stepwise regression procedure 
involves examining at every step the variables incorporated in the error model 
in previous steps. A specific variable is deleted from or entered into the 
model by using the Gaussian Elimination Method for solving the linear system 
of normal equations in the regression. 

Although C-band radar error models are used in the development, the 
procedures can be adapted to other types of tracking systems. Results from 
application to the Apollo-Saturn 501 (AS-501) flight test data are presented and 
indicate generally acceptable truncated C-band radar error models. 


INTRODUCTION 


The errors in a given measurement obtained from a tracking system can 
be classified as random and systematic. The random errors are random in 
character and thus cannot be predicted. The effects of random errors can, 
however, be minimized by using a sufficiently large number of observations in 
a least squares reduction. The combined effects of sources of significant 
systematic error on each quantity measured by a tracking system are frequently 
represented by analytical expressions referred to as error models. The 
problem of evaluating these systematic errors represented by error models is 
of importance in determining an accurate flight trajectory from the basic tracking 
measurements. 



A method for accomplishing this evaluation is provided in TEMS, an 
acronym for Tracking System Error Model Studies. The overall objectives of 
the studies are twofold: 

( 1) To evaluate systematic errors in the tracking system measure- 
ments used in determining a postflight trajectory. 

(2) To conduct analyses to establish truncated tracker error models 
to represent the systematic errors. 

Basically, the TEMS Multiple Regression Analysis Method [1, 2, 3] involves 
establishing the tracker errors and then determining, in the least squares 
sense, error model expressions to describe these established errors. Results 
from application of the method to Apollo-Saturn flight test data are presented in 
Reference 4. In these references, truncating the total error models was 
required because of highly correlated coefficients that made an insignificant 
error contribution. The approach to constructing truncated error models using 
the TEMS results has been based on the significance of an individual variable 
and its correlation with other variables. This approach actually constitutes a 
qualitative examination of a subset of regressions from the all possible 
regressions approach. 

There are other model building approaches. These include: backward 
elimination procedure, forward selection procedure, stepwise regression 
procedure, and stagewise regression procedure. These are not unique, however, 
so far as selecting which of several independent variables should be used in a 
regression equation to describe a response variable. The different procedures 
mentioned above do not necessarily give the same results when applied to the 
same problem. After a careful study of all the procedures, it was concluded that 
the stepwise regression procedure had the greatest potential for useful appli- 
cation to the TEMS error model construction problem. Other applications of 
the stepwise procedure can be found in References 5 and 6. Additional infor- 
mation on the other procedures is also contained in Reference 5. 

The basic approach in the Stepwise Regression Analysis involves 
examining at every step the variables incorporated into the regression model 
in previous steps. At a given step in the analysis, a specific variable is 
deleted from or entered into the regression model by using the Gaussian 
Elimination Method for solving the linear system of normal equations. 

This method is ideally suited for application of the stepwise procedures 
since it obtains elements of the solution vector one at a time. The final 
regression model results in only the most significant variables being 
retained in the model. 
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Application of the TEMS method and the Stepwise Regression Analysis 
to C-band radar tracking systems operating on the Apollo-Saturn V AS-501 
flight is presented later in the report. A summary of the C-band radar 
tracking data utilization on the first burn flight phase (launch to parking orbit 
insertion) and the S-IVB second burn flight phase (S-IVB reignition to S-IVB/ 
CSM separation) is included. The approach given by the guidelines in this 
report for obtaining truncated error models to describe the systematic 
errors has generally resulted in acceptable models for the AS-501 first and 
second burn data. Several cases are noted on the stepwise results for the AS-501 
data where the introduction of additional variables into the regression does not 
improve significantly the cr ■ curve fit value. It appears that a rather critical 

examination of results from application of the stepwise procedures is required 
in order to obtain meaningful and useful information for input to the TEMS 
program . 
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THETEMS MULTIPLE REGRESSION ANALYSIS METHOD 


Introduction 


The development of mathematical procedures in the TEMS Multiple 
Regression Analysis Method is presented in this section. The development 
provides for a comprehensive evaluation of systematic errors in measure- 
ments obtained from various radar tracking systems. Detailed mathematics 
for a rigorous least squares adjustment of the radar error model parameters 
(coefficients) is presented. Included are provisions for the use of a priori 
information about the coefficients of the error model expressions. Coefficient 
variances can also be constrained to be consistent with a priori values. In 
addition, a priori estimates for the coefficients can be entered into the adjust- 
ment. Functional relations between the coefficients are also considered in 
the development. An important by-product of the adjustment is the variance- 
covariance matrix of the estimated error model coefficients. 

The C-band radar tracking system error models used in the development 
are discussed in Appendix A. It should be noted, however, that the development 
of the method for application to tracking systems other than radars is analogous 
to that in this section. A similar development for the AZUSA (Glotrac Station I) 
tracking system is contained in Reference 3. 

Observational Equations for the Least Squares Adjustment 
with Parameter Constraints 


The basic tracking data which serve as input consists of the radar 
measured tracking parameters (R°, A 0 , E°) and of a reference trajectory 
representing the best estimate of the trajectory from a composite of data. 
These composite data are from various tracking systems such as AZUSA, 
Radar, ODOP, Glotrac, and fixed cameras. The reference trajectory is in 
an earth-fixed plumbline coordinate system (X , Y , Z ) with origin at 

6 G C 

the launch site. These data are transformed into radar reference tracking 
parameters by the two transformations: 
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(a) (X , Y Z ) 
e e e 


(X, 


es 


Y es’ 


Z 

es 


) 


(b) (X Y , Z ) -~(R r , A r , E r ) 
es es es 

Detailed mathematics involved in these transformations are given in 
Appendix B. The tracking errors for the particular system under consideration 
and for the i-th observation are thus determined from the equations: 



An additional transformation of the form: 

,C) (X e S ’ Y es’ V ~ (X ’ Y ' Z) 

is required to determine the position of the vehicle in an earth-fixed ephemeris 
coordinate system with origin at the tracking site. This transformation is 
also given in Appendix B. 

The fundamental observational equations are expressed in the following 
form (i = 1, 2 n) : 


AR® - AR. - V_ = 0 ^ 

x l Ri 

AA? - AA. - V A . = 0 > 

l l Ai 

AE° - AE. - V„. = 0 

i l El / 


( 2 ) 


The functional expressions for the systematic errors in equation ( 2) are 
given by (i = 1,2 n) : 
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AR t = C 0 + CjRj + + C 4 (-. 022 cosec E.) + C 5 (X./R.) 

+ C e (Y i /R.) + C 7 ( Z./R.) + Catj 

AA^ = D 0 + D 7 A. + D 3 A i + D 5 tan Ej + D e sec E. + D 7 tan E. sin A. 
+ D 8 tan Ej cos A ^ + D 9 [(sin A^ cos A.)/XJ 
+ D 10 [(-sin A i cos A.J/Yj] + sec E i 
AE a = F 0 + FjE^ + F 3 E + F 5 ( -sin A^ + F e cos A . 


7 l R i si “ E i 


- 10 _e 1 cotanE. + F*[(-X. tanE.)/R?] 


l l l 


+ F 10 [(-Y i tan E.)/R*] + F 11 [(cos E.)/R.] +FJ 2 E.COSE. 



The following constraints in the form of functional relations between the 
coefficients are imposed upon equations ( 3) : 



Equations (3) are thus rewritten as (i- 1,2 n) 
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where 


AE. = C 0 + + C 2 ? 2 . + C 4 f 3 . + C 5 r 4i + C^. + C^.-f C 8 ^. 

AA. = D 0 + Cga^ + D 3 a 2i + D 5 a 3 . + Dga^ + D^a^ + D 8 a 6i 
+ C 5 a ?i + C e a 8i + D 11&9i 

AEj = F 0 + C 2 e li + F 3 e 2i + F > 8 e 3i + D 7 e 4i + ^ e 5i + ^5® 6i 
+ C«e 71 + C,e 81 + Fj2e gl 


esj 


r li = B i 

r 4i 

r 9 . = R. 

2i i 

? 5i 

r n . = 022 cosec E. 

3i l 

r 6i 


/v/ 

r 7i 

a n =A i 

6 ii 

* 2i = A , 

e 2i 

a„. = tan E. 
3i i 

®3i 

a.. = sec E. 
4i i 

e 4i 

a_. = tan E. sin A. 
5i l i 

®5i 

a„_. = tan E. cos A. 

6i i i 

a^. = ( sin A. cos A. ) /Xj 

e 6i 

a g . = (-sin A. cos A^/Yj 

®7i 

• 

a^ = A. sec E. 
9i l l 

e 8i 


®9i 


= X./R. 

1 l 

= VR, 

" Z i /E i 

= 'i 

= E. 

1 

• • 

= E. 

l 

= -sin 
= cos Aj 

= [ • ■ ° 22 — - 10-1 cotan E - 

[_R sin E i J l 

= ( -X. tan Ej) /R? 

= ( -Y. tan E.)/R? 
l i i 

= (cos E.)/R. 

= E. cos E. 

l l 



The coefficients in equation (5) are written as: 


C 0 - Co ^ Cq 
C! = Cj + 6Ci 


F12-FJ2 + 6F12 7 


where C 0 , C ls . . . , F 12 are approximations (which can be chosen equal to zero 
if desired) to C 0 , C ls . . . , F 12 , respectively; and <5C 0 , <5C ls . . . , 6F 12 are 
unknown corrections which will be determined in the least squares sense. The 
systematic errors can then be written as (i =1,2, . . . , n) : 


AR. = AR. + 5 AR. 

li l 

AA. = AA. + <5 AA. 

li l 

AE. = AE. + 6 AE. 

li l 


(7) 


where 


^ r^* ^ ^ 

AR. = C 0 + C 1 r 1 . + C 2 r 2 . + . . . + C 8 r 7j . 
AA. = D 0 + C^. + D 3 a 2 . + . . . + D^. 


AE. - F o + c 2 e li + F 3 e 2 . + . . . + F i 2 e 9i 


( 8 ) 


8 



and 


6AR. = 6C 0 + 5C 1 r 1 . + 6C 2 r 2 . + . . . + SC 8 r 7 . 

<5AA. = 6D 0 + 6C2a . + 5Dga + . • • + SD^a . 

l li 2i 9i 

6AE. = 6F 0 + 6C 2 e 1 . + 6F 3 e 2 . + . . . + SF^. , 


(9) 


Substituting equations (8) and (9) in equations (7) and the result in 
equations (2) yields the following observational equations (1 = 1,2 n) : 




AH. - [ ( C 0 + C-ir^ + . . . + c 8 r 7i ) 


+ (6C 0 + 6C 1 r li + . . . + 6C 8 r 7 .) ] - V R . = 0 


AA° - [ (D 0 + C 2 a 11 + . . . + D^a^) 

+ (SD 0 + <5C 2 a ii + . . . + 6D lt a 9i ) ] - = 0 




AE° - [ ( F o + C 2 e 11 + . . . + F 12 e g .) 


+ ( 6F 0 + 6C 2 e^^ + . . . + <5F 12 e g ^) ] - - 0 


The system given by equations ( 10) can be written in matrix notation as: 


( 10 ) 


N - (BC + B 6) - V = 0 


(ID 


9 



m i 



10 



ir 



11 




12 


( 91 ) 



In addition to the observational equations (11) , error model coefficient 
observational equations are introduced into the adjustment using the following 
equations: 


Cn = Cn +V, 


Ci = Cl + V 


Cl 


F 12 — F }2 + V 


F 12 


(17) 


where the superscript °° denotes an a priori value of the coefficient and 

V , V , . . . , V are the corresponding observational residuals. 
c o Cj F 12 

Substituting equations (17) into equation (6) and rearranging these equations, the 
following observational equations for the error model coefficients are obtained: 


V p - <5C 0 
'-'0 


Cn + C n =0 


- SCi - C t + C t =0 

C 1 


V F - 5Fi2 - FJ2 + F 12 - 0 
r 12 


(18) 


Equations (18) are written in matrix notation as: 


V - 8 - e = 0 

where 


V 

[18x1] 


1 

• < < 
o o 

h* o 

1 

(20) 


CO 

C 0 - C 0 

OO 

Ci-Ci 

V F 
V 12 


e = 

[18x1] 

^ * OO 

F 12 “ Fi2 


(19) 


( 21 ) 
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The overall set of observational equations is then given by: 


*\ 


N - (B C + B 5) - V = 0 
V- 6-7=0 




J . 

The Normal Equations 


(22) 


The final system of observational equations given by the matrix equations 
(22) consists of (3n+ 18) linear equations in the 18 unknown corrections to the 
approximated error model coefficients. In general, for the application herein 
the total number of observational equations (3n+ 18) is large compared to the 
number of unknown corrections. Thus the system consists of an over-deter- 
mined set of equations wherein the effects of random observational errors can 
be minimized in a least squares reduction. According to the principle of least 
squares , the best representation of the data is that which makes the weighted 
sum of the squares of the residuals a minimum. Thus the function to minimize 
is: 


n 


n 


f(6C„, dCj SF„) « 2( W 3 ,-2 V R, )+ 2< W 3i-l V il> 


i=l 


i=i 


Z< W 3i V L )+W > V C 0 +W * V C. + '-- +V 
i=l u 1 


W 18 V 


(23) 
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The condition which fulfills the minimizing requirements is that the 
partial derivatives with respect to each of the 18 unknown corrections be zero. 
Thus the following 18 normal equations are obtained: 


8f 

n 

v 

r av i 

W V R1 

n 

, v 

r av Ai i 

TTT T7 

) 

3(6C 0 ) 

L 

i=l 

3i-2 Ri 3(6C 0 ) 

+ L 

i=l 

3i-l Ai 3(6C 0 ) 



n 

+ 2 
i=l 


W M V. 


9V Eil 


31 Ei 3 ( 6C 0 ) 


9V C 0 

+ w,v - — — - = 0 

WlV C 0 3(6C 0 ) 


> (24) 
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Note: Equation (24) is concluded on following page. 
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0 0 0 0 0 

cP'q/ 0 0 0 0 

0 oVojgi 0 0 0 (26) 

0 0 a 2 / cr^n 0 0 

0 0 0 oVo^n 0 

0 0 0 0 oV a^ n 
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p 

[ 18x3n] 


P = 
[18x18] 


W = 
[ 18x18] 


3V ri 

8v ai 

9V E1 

dV„ 

Rn 

9V a 

An 

aV En 

3(60,) 

9(<5C 0 ) 

O 

o 

to 

CO 

” 9(6C 0 ) 

9(6C 0 ) 

9(5C 0 ) 

8V ri 

CD 

>•* 

h-*- 

aV Ei 

9V„ 

Rn 

9 V . 
An 

aV En 

8(6C,) 

9(6Ci) 

9(6Ci) 

" 9(5Ci) 

9(6Ci) 

9( SC*) 

3V ri 

aV Al 

aV El 

dV„ 

Rn 

9 V . 
An 

9V,-, 

En 

9(6C 2 ) 

9(6C 2 ) 

9(6C 2 ) 

9(6C 2 ) 

9(<5C 2 ) 

9( (5C 2 ) 

• 

aV Rl 

aV Al 

9V E1 

aV T, 

Rn 

9V a 

An 

9V__, 

En 

9(6F 12 ) 9(<5F 12 ) 9(<5Fi2) 

" 9(6F 12 ) 

9(6F 12 ) 9(<5F 12 ) 


(27) 


av r 

u o 

9(5C 0 ) 

0 


0 

aV Cl 

a(<5c7) 


(28) 


ay. 


12 


9(6 F 12 ) 


^<Ar 


c^q/o 2 


0 

0 




12 


(29) 
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It is perhaps important to point out that if the j-th error model coefficient 
is to be unconstrained, then the j-th diagonal element of equation (29) is equated 
to zero. It has also been found from experience that if the error model fits the 
physical situation and there is no significant correlation between the coefficients, 
then erroneous (within reasonable limits) observational weighting in equation (26) 
does not have a significant effect upon the determination of the coefficients. The 
associated estimates of variances and covariances become erroneous, however. 
Additional information on the subject of observational and parameter weighting 
in the least squares adjustment is contained in Appendix C of Reference 2. 

The matrix expressions for the residuals V and V in equation (25) can 
be obtained directly from the observational equations (22) and written as: 


V = N- (BC+BS) ' 

V = 6 + 1 

Substituting equations (30) in equation (25) 
equation for the normal equations: 


(30) 

yields the following matrix 


PW [N - (B C + B <5) ] + PW(e + S) = 0 


(31) 


Solving for <5 : 


6 = (-PW B + PW) -1 [ - PW (N - BC) - P W e] 


(32) 


Detailed examination of B and P reveals that the solution given by equation (32) 
is equivalent to: 

6= (B T WB + W) -1 (B T WN - B T W BC - W e) ( 33 ) 

After 5 has been found from (33) , this value is added to the initial approximation 

rp — 

C to obtain the adjusted error model coefficients. The matrix (B W B + W) 
is the coefficient matrix of the reduced normal equations. The inverse of t his 
matrix multiplied by the unit variance provides an estimate of the output variance-co- 
variance matrix of the error model coefficients resulting from the regression. 
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Properties of the Adjustment 


The vector of least squares residuals for the observations may be 
computed directly from the observational equations or from equations (30) 
after the normal equations have been solved for 6. The least squares 
residuals thus obtained can be rewritten in the more familiar quadratic form 
as: 


s = v T wv + v T wv 


(34) 


An estimate of the unit variance resulting from the adjustment is given by: 


a o = S/f (35) 

where f is the degrees of freedom involved and is the number of observations in 
excess of the minimum required for a unique solution. 

This is given by: 

f = 3n (36) 


where 3n is the total number of observations less the number of unknowns. 

An estimate of the variance-covariance matrix for the parameters re- 
sulting from the regression analysis is given by: 



For the 18 x 18 case, the matrix cr is given by: 

c 


a c = 


at 


C 0 Ci 


a 


c<A 

2 


C 0 F 12 


cr 


'i 


CiF 12 


( 38 ) 


L a a a ! 

C 0 F n CjF 12 F ffl J 

The elements of the covariance matrix of the error model coefficients thus 
per mits the accuracies of the final results to be estimated. 


MATHEMATICAL DEVELOPMENTS IN A STEPWISE 
REGRESSION ANALYSIS 

Basic Stepwise Approach 


The basic approach in a stepwise regression analysis involves examining 
at every step the variables incorporated into the regression model in previous 
steps. A variable determined to be significant at an earlier step may be 
insignificant at a later step because of its relation with other variables in the 
current regression. Results from a given step in the analysis provide statistical 
F tests whereby it can be determined if a specific variable should be deleted from 
or entered into the regression model. The test for deletion is made before a 
variable is added to the regression. The procedure is terminated when no 
more variables will be admitted and no more will be rejected. A summary of 
the basic stepwise approach is given in Figure i. 


19 




FIGURE 1. BASIC STEPWISE APPROACH 

Regression Analysis for Data Centered 
about the Mean 


Assume that the observed response variable Y° is to be estimated by 
the model: 


Y. = b + b Z.. + b Z„. + 
1 0 1 li 2 2i 


+ b Z . 
P Pi 


(39) 


where i = 1, 2, . . . , n. This model is rewritten as: 


Y i - V »i Z li + <b l Z 1 - b l V + b 2 Z 2i + <»2 Z 2 - b 2 Z 2 > <40) 


+ ... + b Z . + (b Z -b Z) 

p pi p p p p 
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Regrouping the terms in equation (40) , we obtain the following equation: 


Y i = <b 0 + b l Z 1 + b 2 z 2 + ••• + b p V + b l (Z li- v 


+ b <Z - ZJ + ...+b (Z.-Z) . 

2 2l 2 p pi p 


( 41 ) 


If we let: 


z,. = Z.. - Z. 
li li 1 


Z 2i Z 2i ” Z 2 


z . — Z . — Z 
Pi Pi p 


b 0= V b l Z l +b 2 V-" tb p 2 p 


(42) 


Then equation (41) can be written as: 


Y. = b' + b z + b z_. + 
i 0 1 li 2 2i 


. . + b z . 
P Pi 


(43) 


At this point the input data matrix of observations appears as shown in 
Table I. 


TABLE I. INPUT DATA CENTERED ABOUT MEAN 
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The sample estimates b^, b^, b , 


b are obtained by minimizing the 


weighted sum of squares of deviations between the observed and predicted 
response values. Since observations of the response are assumed to have 
equal reliability, unit weighting factors are used. The minimizing function is 
then of the form: 


n 

M = 2 (Y° - Y°) 2 

i=l 1 1 


(44) 


The requirement for M a minimum yields the following normal equations: 


8M 

9b 0 


= 0 


9M 

9b. 


= 0 


9M 

9b 


= 0 


(45) 


Note the exact form of the first normal equation obtained by the indicated 
partial differentiation: 


9M 

9b' 

0 


n 

2 (Y° - b' - b z.. - bz 

. , i 0 1 li 2 2i 

i=i 


b z .) 
P Pi 


(- 1 ) = 0 


n 


n 


n 


n 


= - 2 Y7 + nb» + b, 2 z + b 0 2 z + ... + b S z . = 0 
. . l 0 1 . , li 2 2i p . . pi 

1=1 1=1 1=1 1=1 


-nY + nb' + nb z +nb z„ + ... + nb z =0 

0 11 2 2 p p 


Y + bl + b, z, + b n z 0 + . . . + b z =0 

0 112 2 p p 


(46) 
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But: 


Z l = 


n 

2 

i=l Z li 
n 


2 (2 u- z t> 

1=1 

n 


^ Z 11 - Z l^ + ^ Z 12 ” Z l^ + * * ' + ^ Z ln - Z l^ 


n 


n 

= S Z n - nZ 1 = 0 , (47) 

i=l 

n n 


and the same holds for z^, 



z . 
P 


Thus, equation (46) becomes: 


b o 


Y 


(48) 


Substituting equation (48) into equation (42) , the constant term in 
equation (39) is given by: 


b o = Y 


b l Z 1 - b 2 Z 2 


b 

P 


Z . 
P 


A similar substitution in equation (43) yields: 


(49) 


c 

Y. 

l 


Y + 


b z + 
1 li 


b„ z_. 

2 2i 


• • « b z • 
P Pi 


(50) 


This is an alternate form to equation (39) and we have to estimate one less 
parameter since equations (48) and (49) hold regardless of the values for 
b , b , . . . , b in equation (46) . The remaining normal equations are given 

1 ^ P ^ 

by: 
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9M 

8b, 


n 

2 

i=l 


z,.(y. -bz -bz -...-b z.) 
lx J i 1 lx 2 2i P pi 


= 0 


8M 

9b„ 


n 

2 

i=l 


z„. (y. - b, z.. - b„ z„. 
2x '°i 1 lx 2 2x 


b z .) 
P Pi 


= 0 


( 51 ) 


0M 

0b 


n 

2 

i=l 


z . (y. - b, z,. - b z 0 . - . . . - b z .) 
px J i 1 lx 2 2x P pi 


- 0 


Let: 


n 


n 


I = 


S TT = 2 z T . z = 2 (Z -Z)(Z -Z) J 

U . , Ii Ji . , Ii I Ji J 
i=l 1=1 


1 , 2 , 

1 , 2 , 


, p for each 
>P (52) 


n 


n 


S = 2 z y = 2 (Z n - Z ) (Y° - Y) I = 1,2, . . . ,p (53) 
i=i i=l 


s yyL^ < Y ° - Y) < Y i - Y) 

1=1 1=1 


(54) 


Then the normal equations (51) can be written as: 


b l S il + b 2 S 12 + *•* + b p S lp S 1Y 


b l S 21 + b 2 S 22 + * ' * + b p S 2p S 2Y 


(55) 


b, S. + bS+...+b S = S 
1 pi 2 p2 P PP PY 
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The solution of equations (55) for the b's yields: 


1 

o' 

1 


S il 

S 12 ••• 

S lp 

-1 

S 1Y 

b 2 

= 

S 2i 

S 22 •** 

S 2p 


S 2Y 

• 

b 

P 


• 

• 

s < 
Pi 

S p2 

S PP 


♦ 

• 

S 

pY 


(56) 


The estimate of the standard deviation of the response variable is 
obtained from the residual deviations and is given by: 


= 


n 


i=l 


(Y? - Y C ) 2 
1 l 


n - d 


1/2 


(57) 


It can be shown [7] that the standard errors of the estimated partial 
regression coefficients given by equations (56) are given by: 




(58) 


where the c's are elements in the inverse of the S matrix in equations (56) . 
T 

given by: 


The standard error of the term b^ = Y in the regression equation is 


V = <V / ' 

0 


(59) 
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Confidence intervals [ 8] for the estimated coefficients are obtained from: 



(60) 


The specific value to use for 6 can be obtained from a "t" table containing the 
percentiles of the "t" distribution. It is dependent on the number of degrees of 
freedom (n - d) and the particular confidence limits selected. If 95 percent 
confidence limits are selected, then the confidence intervals given by equation 
(60) can be interpreted as follows. If the experiment is repeated with another 
set of data, then another set of limits would be obtained. Continuing the process 
in this manner, then 95 percent of the limits obtained would cover the particular 
coefficient. 


Centering the input data about the mean as shown in Table I reduces the 
absolute size of the numbers entering the computations. Since the centered 
data are required in obtaining the linear correlation matrix, no additional 
computations are required that would not otherwise be available. 


Regression Analysis for Data 
in Correlation Form 


The regression analysis can be transformed into a form which involves 
linear correlations [ 8] . This is accomplished by making the following trans- 
formations on the centered data of Table I. 



r*j 


z 


2i 


z . = 
Pi 


(Y. - Y) / 'S/S^. 

^rV 7 ^ 

(Z 2i- 


(Z . - Z ) / n/S 

pi p pp 


(61) 


26 




The same transformations on the dependent and independent variables in the 
model given by equation ( 50) yields the following model: 


y. 


1 '^ = b l' / ^7 Z U +b 2'^"2i + 


+ b 's/s z . 
P PP Pi 


(62) 


Or: 


where 


y. = a, z,. + a z„. + ... + a z . 
J i 1 li 2 2i p pi 


(63) 


an, = b 


11 


1 1 / S. 


YY 


o' = b 


22 


2 2 /S. 


YY 


(64) 


a = b 


/S 

! pp 


p p J s. 


YY 


The new coefficients o' , o' .... a are to be estimated from the 

12 p 

transformed data y., z z . . . , z ... The minimizing function is: 


n ~ ~c 2 

V= S (y -yV . 

i=l 


(65) 


The normal equations are: 

av 


da. 


= 0 


av 

act” 


= o 


( 66 ) 


av 

act* 


= o 
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Evaluation of these equations yields: 


U V _ ~ ~ ~ ~ 

— — = 2 z A . (y. - a. z . . - O' z„. - . 

3 a. . , ii J i 1 ii 2 2i 

1 1=1 


. . -Q! z .) = 0 

P Pi 


2 z„. (y. - a. z.. - cr z„. - ... - a z .) = 0 

i=1 2i •’l 1 li 2 2i p pi 


— — = 2 I z . (y. - ol. z.. - O' Z-. - . . 

3o! . . pi i 1 li 2 2i 

P 1=1 1 


-a z .) =0 
P Pi 


S = 2 z z I = 1, 2, . . . ,p for each J = 1,2, . . . ,p (68) 
U . , 11 Jl 

1=1 


S IY 2 Z Ii y i ’ I - 1,2, . . . ,p . 
1=1 


Then equations (67) can be written as: 


oi S + oi S + . . . + a S = S,„ 
1 11 2 12 p lp 1Y 


oi S + o' S + . . . + a S„ = S.__ 
1 21 2 22 p 2p 2Y 


a. S . + O' S _ + + a: S = S 

1 pi 2 p2 ”• P PP pY 
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The solution for the a's becomes: 



Examination of equation (68) reveals that: 

g - !a 

M Vs“s 

II JJ 


1=1 

rs. ~ 9 H = 2 ~ 

/ 2 ( z - Z ) 2 (Z -Z) 

\/ . , ' Ii r , v Ji j' 
1=1 1=1 




= r 


IY • 


( 71 ) 


( 72 ) 


( 73 ) 
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Thus, the S elements of equations (71) are the linear correlation coefficients 
for the independent variables and the response. Equations (71) can then be 
written as: 



(74) 


The constant term is obtained by solving equations (64) for b^, b^ b 

and substituting in equation ( 49) : 



The estimate of the standard deviation of the transformed response 
variable is given by: 


5, ~ c 2 

2 (v. - y ) 

i=l 

n - d 


The standard errors of the new regression coefficients are given by: 
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where the c's axe elements In the inverse of the r matrix in equation (74) . 
The standard error of the term = Y is obtained by first writing equation 

(85) , which appears later in the discussion, as: 

a Y = ^ S YY °Y • (78) 


Using this result in equation (59) , we see that: 

<r h , = °Y ^YY (79) 

0 \/n 


This equation is necessary if the input data are in correlation form. Equation 
(59) , however, is sufficient when the input data are centered about the mean. 


Alternate Form for the Standard Deviations of the 
Estimated Partial Regression Coefficients 

The standard errors of the estimated partial regression coefficients 
as given by equations (58) can be expressed in another form. Equation (76) 
is first written as: 


~2 


o 


Y 


~ ~c,2 ~c.2 ~ ~c.2 

( Yj -Yj) + (y 2 -y 2 > + ••• + <y n -y n > 

n - d 


(80) 


Substituting z , z^, . 


z . of equations (61) in (63) and the result in 
Pi 


equation (80) and y^ of equations (61) in equation (80) we see that: 
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L 



~2 __ 1 
°Y ~ n-d 


(Y° - Y) 


-> 2 


“i (z ii- z i> 


Qi ( Z - ZJ 
2 ' 21 2 




YY 


■sfs 


11 




22 


y v - v 


'(y° 2 -y) 


a i< Z «- 2 l> 


“«(Z„-ZJ 
2 22 2 


a (Z -Z ) 
P P2 P 


va 


YY 


n/s" 


11 


\Ts 


22 


«/s" 


PP 


+ . 


(Y°-Y) 


“l^ln- 2 !* 


's/s" 


YY 


N/S 


11 


y Z 2 n- Z 2 > 


Q! (Z - Z ) 

■- P . P n E_ 


-Vs" 


PP - 1 


(81) 

Substituting equation (64) for the standard partial regression coefficients in 
equation (81) and rearranging the equation, we find: 


_2 


-2 111 


cr = 


Y n-d S 


YY 


^-Y)-b 1 ( Zll -Z I )-b ! IZ ai -Z 2 ). 


-b <Z -Z ) 
p pi p 


YY 


-.2 


< ^2- 5) - b 1 (Z 1 2- 2 l ) - b 2 (Z 22- Z 2 ) --”-y Z p2-y 


+ . . . + 


YY L 


(Y n" Y) " b l (Z m- V-V Z 2n- V 


-b ( Z - Z ) 
p pn p 


( 82 ) 
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c — 

Then, by substituting (Y. - Y) as obtained from equation (50) in equation 
(82) , the following equation is obtained: 


~ 2 
cr = 


Y (n-d) S 


YY L 


(Y° - Y) - (Y° - Y) 


(Y° - Y) - (Y° - Y) 


+ . . . + 


(Y° - Y) - (Y° - Y) 
n n 


n 


c ,2 


(n-d) S. 


YY 1 i=l 


2 (Y~ - Yp 


(83) 


But from equation (57) , equation (83) can be written as: 


V s 


YY 


(84) 


Or: 


cr = 


o* / n/s 
Y YY 


(85) 


This equation gives the relation between the response variable and the transformed 
response variable. 

This result can be substituted in equation (77) for a typical j-th 
coefficient: 


cr 

oi. 

1 




(86) 


Since a. is a linear function of b., then the variance in a. can be written as: 

i A ' A 


j 


j 


j 



(87) 
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By substituting <r from equation ( 87) in equation ( 86) , the following 

j 

equation is obtained for the standard error of a typical j-th partial regression 
coefficient: 


r. = , Y ■ n/cC . (88) 

b. n/s H 

J YY 


Equation ( 88) is used in place of equation ( 58) when the regression analysis 
has been transformed into the form involving correlations ( input data in 
correlation form) wherein the inverse elements of the S matrix are not 
available. 


Alternate Form for the Standard Deviation 
of the Response Variable 


The residual sum of squares can be written as (see Reference 5) : 


n 

2 

i=l 


( Y° - Y?) 2 
1 1 


n 

2 

i=i 


(Y°- Y) 2 
x 


n 

2 

i=l 


( Y? - Y) 2 . 

l 


(89) 


That is, the residual sum of squares is equal to the total sum of squares about 
the mean minus the sum of squares due to regression. But from the definition 
of the multiple correlation coefficient: 


n n o n o—2 

2 (Y_ - Y7 = 2 ( Y, - Y) b: 


Y. 12... p 


i=l i=l 

Substituting equation (90) in equation (89) , we find: 


(90) 
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n 

2 

i=l 


(Y°- Y? ) 2 
1 l 


n 

S 

i=l 


(Y.°-Y) 2 



12 ... 


P 


YY 




Then by substituting equation (91) in equation (57): 


= 


YY 


4 

n - d 


(l - R". , 0 ) 

' Y . 12. . .p / 


i 

2 


( 91 ) 


(92) 


This expression enables the estimate of the standard deviation of the response 
variable to be expressed in a more convenient computational form. Its use has 
an advantage over equation (57) when the Gaussian Elimination Procedure is 
used to solve the normal equations. This is discussed later on in this 
report. 


Partial Correlation Coefficients 


The partial correlation coefficient of the variable Z (not in the regres- 

Si 

sion) with the response Y — given that Z , Z , . . . , Z are already in the 

12 p 

regression — is denoted as p „ . This statistic is used in the step- 

qY . 12. . . p 

wise procedure to determine which of several variables not in the regression 
should be considered for entry into the model at a given step. The variable 
considered for entry is the one whose partial correlation with the response is 
highest. 

Mathematically, p is found by determining the simple linear 

correlation between the residuals from the regression: 


Y. = b + b Z + b Z + ... + b Z. 
l 0 1 li 2 2i p pi 

and the residuals from the regression: 


(93) 
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z c . = f>. + 6, z, . + £ 0 z 0 . + . . . + 6 z . 

qi 0 1 li 2 2i p pi 


( 94 ) 


Let the residuals from these regressions be denoted as: 


o c 
VR, . = Y. - y: 
li i i 


VR n . = Z . - Z . 
2i qi qi 


Then: 


(95) 


qY . 12. . . p 


2 <VB a VB 21 > 

1=1 

/n 2 n 2 ~ 

/S VR,. S VR* 

/. . li . , 2i 

/i=l i=l 


(96) 


An alternate form for the partial correlation coefficient involves the 
linear correlation coefficients. Let: 


— 




— 1 

-1 

— 





r ll 

r !2 •" 

r ip 

r iq 

r iY 


C 11 

C 12 

C lp 

°iq 

C 1Y 

r 

r ... 

r 

r 

r 


c 

c ... 

c 

c 

c 

21 

22 

2p 

2q 

2Y 


21 

22 

2p 

2q 

2Y 


• 

• 

• 

• 
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• 


• 


• 

• 

• 

• 

• 


• 

• 

• 

• 

• 

• 

# 

• 

# 

• 


• 


• 

• 

• 

r . 

r „ ... 

r 

r 

r 

— 

c , 

c ... 

c 

c 

c __ 

pi 

p2 

PP 

pq 

py 


Pi 

P2 

PP 

pq 

pY 

r , 

r* 

r 

r 

r 


c . 

c _ ... 

c 

c 

c 

qi 

q2 

qp 

qq 

qY 


qi 

q2 

qp 

qq 

qY 

r 

r ... 

r 

r 

r 


c , 

c ... 

c 

c 

c 

Y1 

Y2 

Yp 

Yq 

YY 


Y1 

Y2 

Yp 

Yq 

YY 


( 97 ) 
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Then from Reference 9: 


p qY . 12. . . p 


-c 


SL Y. 


VC c 

qq YY 


(98) 


Other partial correlation coefficients can be determined having once 
determined the inverse equation ( 97) . For example: 


P 1Y. 23. . . qp 


-c 


1Y 


^ C u °YY 


(99) 


Significance of the Estimated Regression Equation 


The significance associated with an estimated regression equation can 
be considered from an analysis of variance viewpoint. If the estimated 
regression equation is of the form given by equation (50) , then the analysis 
of variance is summarized in Table II. 


TABLE II. ANALYSIS OF VARIANCE 


d f 

Type 

Variation 

SS 

MS 

F Value 

n - i 

Total 

n o-2 

S = S (Y -Y) 

i=l 



n - d 

Residual 

S(RES) = S ( Y?- Y?) 2 

i=l 1 1 

S(RES) 


1VI\ X\ JliO ) — 

n - d 


d - 1 

Regression 

S(REG) = 2 (Y?- Y) 2 
i=l 1 

M(REG) = S( ^ EG ~ 
d-1 

M(REG) 

M(RES) 


df = degrees of freedom 

SS = sum of squares 

MS = mean square = SS/ df 

n = number of observations 

p = number of independent variables 

d = p+ 1 
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From equation (89) , the total sum of squares can be written as: 


n 

2 

i=l 


(Y°- Y) 2 
1 


n 

2 

i=l 



Y i )2 + 


n 

2 

i=l 



= S(RES)+ S(REG) 


( 100 ) 


That is, the total sum of squares about the mean is equal to the sum of squares 
about regression (residual) plus the sum of squares due to regression. This 
shows that a way of assessing how useful the regression equation will be for 
prediction purposes is to see how much of is given by S(RES) and how 

much is given by S(REG). If S(RES) = 0, then the actual observations are 
described exactly by the regression equation (50) . The ratio defined by: 

4.U...P = s < reg)/s yy < 101 > 


is known as the multiple correlation coefficient and measures the closeness with 
which the regression equation describes the observed data. 


Another useful ratio in regression analysis is given by the following 
equation (Reference 5): 

S(REG) / (d-1 ) 

S(RES) /(n-d) * 


F = 


( 102 ) 


It follows an F distribution with (d-1) and (n-d) degrees of freedom. This 
quantity is used to determine the statistical significance of the regression 
equation under consideration by comparing it with the appropriate F table 
value. For example, determine the 95-percent percentage points F[ (d-1) , 
(n-d) , 0. 95] in the appropriate F table. If the computed F value given by 
equation ( 102) is greater than the table value, then the regression equation 
is statistically significant. An alternate expression for equation ( 102) can be 
obtained by substituting the residual sum of squares as given by equation (91) 
into equation ( 100) to obtain: 

S(BEG) = SyY - SYY (l-R|_i 2 ___ p ) . U03) 
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Substituting equation (91) and equation (103) into equation ( 102) , we see that: 


F = 


Y. 12. 


(n-d) 


(1 - R Y.12...p> <d ' 1) 


(104) 


This equation is a more convenient computational form for the F quantity when 
performing a stepwise regression analysis. This is because the quantity 
2 

( 1 - R ) is available at a given step in the analysis. This is discussed 

X. IZm • *P 

later on in this report. 

It can be shown (Reference 10, Chapter 14), that the F distribution with 
1 and (n-d) degrees of freedom is equivalent to the t 2 distribution with (n-d) 
degrees of freedom. That is, 


F [ 1, (n-d) ] = t 2 [n-d] . (105) 

The significance of the individual partial regression coefficients (b ,b , . . . ,b ) 

1 4 P 

is determined from ^Reference 10, Chapter 2(|t 

\ = V CT b (106) 

i 

which is distributed as t with (n-d) degrees of freedom. By virtue of 
equation ( 105) , this can be written as: 


F i<OUT> = b f /<T b. < 107 > 

1 

with 1 and (n-d) degrees of freedom. The F value is denoted as F.^ to 

indicate its use in the stepwise procedure for leaving an insignificant variable 
out of the regression equation. If the computed F j(q UT j value given by 

equation (107) is greater than the appropriate table value, then the i-th 
coefficient is statistically significant. If the regression analysis has been 
transformed into the form involving correlations (input data in correlation 
form) , then the appropriate form of equation ( 107) to use is obtained by 
substituting equation (92) in equation (88) to obtain for cr^ : 
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(108) 



S YY ( 1 R Y. 12. 
n - d 


i 

~l 2 



The partial regression coefficients can be expressed in terms of the stand- 
ard partial regression coefficients a. from equations (64). This result and 
equation (108) are substituted in equation (107) to obtain: 


F i(OUT) 



(n - d) 


c..(l - R 

u 


2 

Y. 12. 



(109) 


This is the form of equation (107) to use when the input data are in correlation 
form. 


It was pointed out earlier in this report that of several variables not 

in the regression, the variable considered for entry is the one whose partial 

correlation with the response is highest. If this is determined to be 

p (variable Z not in the regression) , then the significance of 

q Y . 12. . . p q 


P qY. 12... p 


is determined from (Reference 10): 


t = 


p qY. 12. . . p 




2 

P qY. 12. . . p 


( 110 ) 


which is distributed as t with (n-d) degrees of freedom and where d is the 
total number of variables. From equation ( 105) this can be written as 


F q(IN) 


2 

p qY. 12. ..p 


(n-d) 


1 


2 

P qY. 12...p 


( 111 ) 


with 1 and (n-d) degrees of freedom. The F value is denoted as F q ( IN ) to 

indicate its use in the stepwise procedure for entering a significant variable 
into the regression equation. If the computed F^^ value given by equation 

(111) is greater than the appropriate table value, then the coefficient corre- 
sponding to variable Z is statistically significant. In the stepwise process 
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the coefficient corresponding to the smallest of qut) * s tested for deletion 

before a new variable is considered for entry. The coefficient corresponding 
to the largest of F is tested for entry after the deletion test. 

The Gaussian Elimination Method as Applied to a 
Stepwise Regression Analysis 


The normal equations given by equations (70) are a system of 
simultaneous linear algebraic equations in a , a , , a . These equations 

1 Lt P 


can be written as: 


a, r, , + a- r, + . .. + a r - r 
1 11 2 12 p lp 1Y 


“l r 2l + “2 r 22 + *'• + “p r 2p r 2Y 


( 112 ) 


a. r. + a'r r + ... + o! r -r 
1 pi 2 p2 P PP PY 


Or in matrix notation: 


r 01 — v 

PP P PY 


where 



“l 



r iY 

a = 

P 

01 2 

(114) 

K! 

11 

r 2Y 


a 

P_ 



V_ 


(113) 


(115) 
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r 


PP 


r ii 

r l2 

r ip 

r . 

r ... 

r 

21 

22 

• 

2p 

• 

r . 

• 

r „ ... 

• 

r 

Pi 

P2 

PP 


(116) 


One method of solving the system given by equations (112) involves 

linear transformations and is referred to as the Gaussian Elimination Method 

[11, 12] . The inverse of the linear correlation matrix r and the solution 
_ PP 

vector a can be obtained from the method. One by-product as applied to the 

regression analysis is the availability of partial correlation coefficients for 

the variables not in the regression. Another advantage involves the estimate 

of the standard deviation of the response variable. Normally this estimate 

is computed from the summation given by equation (57) . A more convenient 

computational form is given by equation (92) . The quantity ( 1 - R^. ) 

Y • • . p 

in equation ( 92) for the p variables in the regression is easily obtained when 
the Gaussian Elimination Method is used. 

Basically the underlying principle of the method involves two elementary 
operations: 

( 1) multiplication (or division) of any equation (any row) by a 
nonzero scalar. 

(2) addition of one equation (one row) to another. 

When these operations are performed on the system of equations (112) , a system 
is obtained which is equivalent to the original system. The basic approach 
is illustrated by augmenting the correlation matrix f as follows: 



— 

PP — 1 


r il r i2 "■ r ip r iY 

1 0 ... 0 

r YY = 

r 21 r 22 r 2p r 2Y 

• * • • 

0 1 ... 0 
• • • 


• * * # 

• * • * 

r r ... r r 

pi p2 pp pY 

• • • 

• • • 

0 0 ... 1 


r t v r 

Y1 Y2 •** Yp YY 

0 0 ... 0 


(p+ 1) x (p+ 1) 

(p+ 1) x p 
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By performing the operations (1) and (2) on rows of this matrix, the 
following arrangement can eventually be obtained: 



(pH) x (pH) (p+.l) x p 

(li8) 


It is seen that appropriate row operations on the r matrix to the left of the 

PP 

line in equation ( 117) produces the identity matrix. The same row operations 
performed on the identity matrix to the right of the line in equation (117) 
produces the inverse of the r matrix. The same row operations on the 

appropriate elements of the (p+ 1) -th column to the left of t^e line in equation 
(117) produces the solution vector a and the quantity ( 1 - R ) . 

P X • LCim • , p 

Elements of the solution vector are actually obtained one at a time. 

This amounts to entering the variables into the regression equation one at a 
time. Thus, the Gaussian Elimination Method is ideally suited for a stepwise 
regression analysis where all p of the independent variables may not be 
desired in the final regression equation. At a given step in the method (before 
equation (118) is obtained) , the array formed by the first p rows and columns 
contains nonzero off-diagonal elements. These nonzero elements indicate 
the variables not in the regression. The partial correlation coefficients of the 
variables, not in regression, with the response can be easily obtained using 
the appropriate diagonal element in this p x p matrix, the diagonal element 
in the ( p+ 1) x (p+ 1) position, and the appropriate element in the (p+ 1) -th 
column. This is clarified in the example in Appendix C. It is pointed out to 
indicate the availability of the partial correlation coefficient required in 
equation (111) for determining This equation provides the test for 

determining the variable to next enter the regression and is discussed further 
in the second example. 
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The following simple example illustrates the basic principles of the 
method. Assume the normal equations are given by: 


r „ 

r 


a 


r 

11 

14 


1 


1Y 

r _ 

r 


a 


r 

.41 

44 


4_ 


4Y 


(119) 


The augmented matrix similar to equation (117) is: 


r 33 


1 

hi 

r i4 

r iY 

1 

0 “ 

r 41 

r 44 

r 4Y 

0 

1 

_ r Yl 

r Y4 

r YY 

0 

0 


(120) 


The nondiagonal elements of the first and second columns must be zeroed out. 
This is accomplished by the following series of operations. Divide the first row 
by r ^, multiply the resulting row by (-r ) , and add the result to the second 


row ( zeros second element of first column) : 

1 r'.. = r../r.. r'.__=r.__/r 


r ' = 
33 


r ' = r /r 
14 14' 11 

r' 

1Y 

r 41 r i4 

-*» t 

1 " " 44 ' r u 

r 

4Y 


11 

r r 
41 1Y 


'll 


Yi 


Y4 


YY 


c ' = l/r 
11 11 


c ' = -r / r 

41 4f 11 




0 



( 121 ) 


Multiply the first row of equation (121) by ( -r ) and add the resulting row to 
the third row (zeros third element of first column) : 
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( 122 ) 


Divide second row of equation (,122) by , multiply the resulting row 
by (-r^ ) > a^d add the result to the first row (zeros first element of second 
column) : 



(123) 

Multiply the second row of equation (123) by (-r.^) add the result to 
the third row (zeros third element of second column): 



(124) 
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It can be verified that equation ( 124) is equivalent to: 



“i 

0 

°T 

no 

C 11 

(NJ — 

°14 

-v _ 

r 33 

0 

1 

*4 

°41 

c44 


0 

0 


'“I 

-a 

4J 


where: 


(125) 


— 


4 

r— 

“ 

r il 

r i4 

”1 

c n 

°14 


r 44_ 


no 

_ C 41 

nsJ 

°44_ 


(126) 


It should be pointed out that the method is reversible. That is, the 
first and third elements of the fifth column in equation ( 124) can be zeroed 
out to obtain equation ( 122). This is especially useful for variable deletion 
purposes when it is desired to leave an insignificant variable out of the 
regression equation. The variable deletion is accomplished by zeroing the 
appropriate elements in the column to the right of the line. The column 
selected corresponds to the particular variable being deleted. This is 
discussed further in the example in Appendix C. 

It is thus seen that zeroing a column in the r matrix results in a 
regression equation with one more variable. Similarly, zeroing a column 
in the c matrix results in a regression equation with one less variable. At a 
given step, the last column of the r matrix contains the regression coefficients 
for the particular variables in the regression. The c matrix contains the in- 
verse of that part of the r matrix corresponding to the variables in the 
regression at a given step. Specific details in the development of the complete 
stepwise procedure are given in Appendix C. 
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Computational Algorithm in the General 
Stepwise Application 


Variable Entry Algorithm . Succeeding elements of the r matrix for a 
variable being entered or deleted from the regression can be generated from 
general algorithms. Let r!. and cl. denote elements of the current matrix. 

If K denotes the independent variable being entered, then the new r!! elements 
can be obtained from the following equations: ^ 


r &5 = r k/ r kK = 1 ; 1 = 1 - K / 

r kj * r k/ r kK '• 1 - K ' J * K 


r n = r T - ( r T r t )/ r 1 

ij ij K iK Kj KK ’ 


elements of 
K-th row 

all other 
elements 


( 127 ) 


The new c! ! elements can be obtained from: 
ij 


°k K = S k K /r kic 1 = ” K 
S kj = J kj /r kK '• 1 = K - i * K 
= (r ; K s ki )/r kK ; 


elements of 
K-th row 


all other 
elements 


( 128 ) 


Variable Deletion Algorithm. An independent variable can be deleted 
from the regression using somewhat similar algorithms. Assume the current 

matrix consists of the r!! and c!! elements. If the K-th independent variable is 

U „ U 

to be deleted, then the new r!!' elements can be obtained from: 
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(129) 


r kk = r ^KK : 1 “ 1 “ K 


ITT — / n T t 


r": = r 


/c" ; i = K, j * K 


Kj Kj KK 
rll' = r!! - (c” r” )/ 


ij 


ij 


iK Kj KK 


elements of 
K-th row 


all other 
elements 


The new c"' elements can be obtained from: 
1J 




c"' 

ij 


- 1 “ - 1 - K ) 
S k/ S KK^= K -^ K ’ 


c.'! 

i] 


- (c! 


.Vo l 


iK Kj KK 


elements of 
K-th row 

all other 
elements 


(130) 


The algorithms given by equations (127) through (130) result from a generali- 
zation of the equations in the Gaussian Elimination Procedure used for solving 
the normal equations. 


Compact Computational Algorithm for Variable Entry and Deletion . 
The order of the augmented matrix at a given step is (p+ 1) x (2p+ 1) . It is 
noted, however, that p of the columns contain zeros and the unity element. 
For every column zeroed in the r matrix, the corresponding column in the 
c matrix contains nonzero elements. The nonzero columns of the c matrix 
can be placed in the r positions corresponding to the column in the r matrix 
that is zeroed out. If this is done from the first step on, then an r matrix is 
set up that can be computed using a compact computational algorithm since 
the other elements are zero and unity and can be generated as needed. For 
example, in equation C-2 of Appendix C, the elements in the ell column 

would be placed in the corresponding positions in the r£ column. The i-th 

element c' would then be designated as r!^. The complete r' matrix for the 

variable in the regression would then be given by: 
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r' 


r' 

ll 

r» 

i2 

r.'o 

i3 

r! 

i4 

r* 

iY 

r H -r i4 r 41 

r i2~ r i4 r 42 

r 13“ r i 4 r 43 

-r /r 
14 ! 44 

r iY _r i4 r 4Y 

r 44 

r 44 

r 44 


r 44 

r 2l" r 24 r 41 

r 22 _r 24 r 42 

r 23 _r 24 r 43 

-r /r 
24 44 

r 2Y _r 24 r 4Y 

r 44 

r 44 

r 44 


r 44 

r 3r r 34 r 41 

r 32 _r 34 r 42 

r 33 _r 34 r 43 

-r /r 
34 44 

r 3Y" r 34 I '4Y 

r 44 

r 44 

r 44 

r 44 

r /r 
41 44 

r /r 
42 44 

r /r 
43 44 

l/r 

44 

r 4Y /r 44 

r Yl" r Y4 r 41 

r Y2 _r Y4 r 42 

r Y3 _r Y4 r 43 

-r /r 
Y4 44 

r YY _r Y4 r 4Y 

r 44 

r 44 

r 44 


r 44 








(131) 

The r!^ elements of this matrix can be computed from the following algorithms 
with K = 4: 


r KK = l/l kK :i= 'i = k 


r Kj r K/ r KK ; 1 K> * * K 


( intersection of row K and 
column K) 

( other elements of row K) / ( 132) 


rig. - -r /r ; i + K, j = K (other elements of column K)| 


r -j - r jj " ( r iK r Ki )/r KK ; 1 * K ’ ^ (all others) 
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At a given step a succeeding matrix (that is, one obtained by entering or 
deleting a variable) can be generated using expressions analogous to equations 
( 132) . The subscript K corresponds to the variable being entered or deleted. 
Note, however, that from the first step on, the c elements must be designated 
as r elements in the manner discussed in the previous example. The algorithm 
given by equation ( 132) for entering or deleting a variable gives the same 
results as equations (127) and (128) for variable entry and equations (129) 
and ( 130) for variable deletion. The obvious advantage in using equations 
(132) is that it involves handling a (p+1) x (p+1) matrix rather than a (p+1) 
x (2p+l) augmented matrix. 


RESULTS FROM THE APOLLO- SATURN AS-501 
VEHICLE FLIGHT TEST 

Introduction 


The Apollo-Saturn AS-501 vehicle was launched at 0700:01 Eastern 
Standard Time on November 9, 1967, from KSC Launch Complex 39, Pad A. 
This section presents results from application of the TEMS Method and the 
stepwise regression procedure to tracking data obtained from this launch. 
Tracking data from seven C-band radars providing coverage on the first burn 
flight phase (launch to parking orbit insertion) and four providing coverage on 
the second burn flight phase (S-IVB reignition to S-IVB/CSM separation) were 
used in th» reduction. The postflight reference trajectory used as the standard 
is presented in Reference 13. 

The relation between the vehicle trajectory for the first burn flight 
phase and the various C-band radar tracking sites is shown in Figure 2. The 
overall first burn tracking data utilization and significant event times are also 
shown in this figure. The second burn flight phase summary with the overall 
tracking data utilization and significant event times is shown in Figure 3. The 
specific tracking data utilization for the first and second burn flight phases is 
shown in Figures 4 and 5. These two figures show the individual time spans 
of usable data as determined from an edit pass through the TEMS program. 
Location data for the launch site and the various tracking stations are given in 
Table III. The minimum elevation angles corresponding to tracking coverage 
at the begin and end times are given in Table IV. 
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FIGURE 4. TEMS AS-501 FIRST BURN RADAR TRACKING 
DATA UTILIZATION 
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FIGURE 5. TEMS AS-501 SECOND BURN RADAR TRACKING 

DATA UTILIZATION 


TABLE IIL LOCATION OF LAUNCH SITE AND C-BAND TRACKING 
RADARS USED IN TEMS AS-501 REDUCTION 



Latitude , 

Longitude , 

Heigh 

Site 

deg 

deg 

m 

Launch Complex - 39, 
Pad A 

28. 608422 

80.604133 

111.58 

Patrick Radar ( 0. 18) 
Merritt Island 

28. 226553 

80.599293 

15.51 

Radar ( 19. 18) 
Grand Bahama 

28.424862 

80. 664404 

12. 02 

Radar (3. 18) 
Grand Turk 

26. 636350 

78.267708 

12. 05 

Radar (7. 18) 

21.462890 

71. 132114 

28.45 

67. 16 (FPS-16) 

32. 348103 

64. 653801 

24. 31 

67. 18 ( FPQ-6) 

32. 347964 

64.653742 

25. 51 

Cape Kennedy (1. 16) 
Antigua Radar 

28.481766 

80.576515 

14.40 

(91. 18) 

17. 144032 

61.792859 

49.37 


a. Elevation above the Fischer Ellipsoid 

b. Elevation of the C-band radar antenna above the Fischer Ellipsoid 
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TABLE IV. ELEVATION ANGLES AT BEGIN AND END TIMES 
FOR FIRST AND SECOND BURN DATA ON AS-501 


Radar 

First Burn 

Second Burn 

Time, sec 

E° , deg 

Time, sec 

E° , deg 

0. 18 

35. 

2. 53 




407. 

7.38 



19. 18 

22. 

2. 07 

11338. 

2.61 


257. 

19. 63 

11840. 

2.85 

3. 18 

99. 

2. 89 

11364. 

1.31 


506. 

2.75 

117.05. 

10.07 

7. 18 

215. 

0. 6 




567. 

0. 9 



67. 16 

379. 

8. 29 




675. 

7. 25 



67. 18 

316. 

3. 87 

11521. 

1.61 


614. 

18. 25 

12386. 

4. 18 

1. 16 

26. 

4. 27 




309. 

15. 29 



91. 18 



11684. 

5. 69 




12380. 

9.26 


Detailed information on the IBM 7094 computer program for application 
of the TEMS method is given in Appendix D. The computer program for 
application of the stepwise regression analysis procedures developed earlier 
in this report is discussed in Appendix E. The utilization of these two pro- 
grams on the AS-501 first and second burn data is illustrated in Figure 6. 
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Stepwise Regression (STEPRG) 
Analysis using results from 
previous block 


Interpretation and analysis of 
stepwise results to determine 
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model 


IBM 7094 


TEMS final pass (or passes) 


FIGURE 6. UTILIZATION OF THE TEMS AND STEPRG 
COMPUTER PROGRAMS 

It should be pointed out that a certain amount of caution has been used 
in the interpretation and analysis of the stepwise regression results. The 
stepwise regression results in Appendix F indicate several cases where the 
0 ^ curve fit value is not improved significantly by the introduction of additional 

variables into the regression. The F and F q UT levels selected for variable 

entry and deletion are also critical in the number of variables entered into the 
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regression model. A value of 3. 5 for and F q UT was arbitrarily selected 

and used in obtaining the stepwise results in Appendix C. A higher value 
would reduce the number of steps and result in a final regression model with 
a lower number of variables. A rather critical examination of the stepwise 
regression analysis results thus seems to be required in order to obtain 
meaningful and useful information for input to the TEMS program. 


Truncated Error Model Regression Analysis Results 


The general approach for obtaining truncated error models to describe 
the AS-501 range, azimuth, and elevation response variables is summarized 
in the following guidelines: 

( i) It was assumed that the survey terms, the rate bias term, and the 
azimuth and elevation velocity lag terms, were not essential in 
obtaining truncated error models to describe the response variables. 

( 2) The first two variables entered in the stepwise regression 
(excluding those left out under the above assumption) were 
selected for consideration in the final TEMS error model. 

( 3) A third variable was considered if an adequate description of 
the response variable was not obtained with the first two or if 
a constraining condition required an additional variable in the 
model. 

This approach actually results in entering the most significant variables into 
the error model. It should be pointed out that the third variable selected 
in guideline (3) often involved selecting one of two variables that represented 
borderline cases so far as the order of entry in the stepwise regression was 
concerned; i. e. , the two variables had p 2 values nearly equal. 

The AS-501 first and second burn truncated error model results obtained 
using guidelines (1) - (3) are presented in Tables V through VIII. Plots of 
the observed and computed response variables and the least squares residuals 
for the truncated models are given in Appendix F. This appendix also includes 
a summary of the stepwise regression analysis results. Coefficient correlations 
from the TEMS program for the truncated models are also presented. The 
average standard deviations of the least squares residuals in Tables V and VII 
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TABLE V. TRUNCATED RADAR ERROR MODEL MULTIPLE REGRESSION 
RESULTS FOR FIRST BURN DATA ON AS-501 VEHICLE FLIGHT TEST 




TABLE VI. COEFFICIENT STANDARD DEVIATIONS FOR TRUNCATED RADAR ERROR 
MODELS FOR FIRST BURN DATA ON AS-501 VEHICLE FLIGHT TEST 


Radar <j For Indicated Coefficient 

K 



C 0 

Ci 

C 2 

C 4 

Do 

D 3 

D 5 

d 7 

Do 

Fo 

Fo 

Terms 

0. 18 
19.18 
3. 18 

0.84 

0.72 

0.36 


0.25E-3 
0.33E-3 
0. 11E-3 

5.21 
3.75 
2. 27 

0. 54E-3 
0. 0010 

0.074 


0.0015 

0.0016 


0.93E-3 

0.0010 

0.24E-3 

0. 100 
0.050 
0. 173 



0. 0024 


y 


0.23E-3 

0. 086 


7 

7. 18 
67. 16 
67.18 
1.16 

0. 95 

1. 15 
0. 84 
1.03 


0. 30E-3 
0.09E-3 
0. 10E-3 
0.89E-3 

1. 97 

0.64E-3 

0.56E-3 

0. 46E-3 

1. 14E-3 

0.971 
0.005 ' 

0.004 
0.073 




0.61E-3 
0. 54E-3 
0.48E-3 
1.57E-3 

6 

0.21E-5 
0. 15E-5 
0.67E-5 



0. 48E-3 

0.014 

0.012 

8 

Q 




0.43E-3 



0. 0025 

O 

n 





l 

No. 













Occur- 

rences 

7 

3 

7 

4 

7 

6 

1 

3 

2 

7 

5 




OS 

o 


TABLE VII. TRUNCATED RADAR ERROR MODEL MULTIPLE REGRESSION RESULTS 
FOR SECOND BURN DATA ONAS-501 VEHICLE FLIGHT TEST 


Radar 






Coefficient 






No. o. 











o 

(J ir A 













VR 

VA 

VE 

Data 


C 0 

c, 


Do 

D 3 

d 7 

D 8 

F 0 

f 3 

m 

Deg 

Deg 

Point:- 

19. 1* 

-25. 89 

-73E-5 


-0.86E-3 

-0. 0117 

-0. 0223 

— 

0. 0018. 

0.4277 

7.30 

0.0050 


492 

; 3. is 

8.46 




0.2989 


— 

-0. 0182 

2. 5233 

7. 16 

0.0038 

0.0055 

322 

1 91.18 

i 

34. 84 

1 -2.43E-5 

: 0.0043 


! 0.4009 


1 — 

| -0.7E-4 

1 -4.1971 

2.22 

0.0038 

0. 0059 

1 684 

67. 18 

8. 83 

-2. 87E-5 


0. 0055 



0. 0012 

-0.0116 


7.28 

0. 0053 

0. 0058 

864 

Average cr 

5. 99 

0. 0045 

0.0061 


TABLE VIH. COEFFICIENT STANDARD DEVIATIONS FOR TRUNCATED RADAR ERROR 
MODELS FOR SECOND BURN DATA ON AS-501 VEHICLE FLIGHT TEST 


Radar 

a.. For Indicated Coefficient 
k 



Co 

c, 

C 2 


»3 

d 7 

o 8 

F 0 

F 3 

Terms 

19. IS 

0. 60 

0.05E-3 

0. 32E--1 

0.33E-3 

0. 057 

0. 0012 

— 

0. 58E-3 

0. 157 

x 

3. IS 

0. 29 


0. a 5 E - 4 

0. 35E-3 

0. 075 

— 

- 

0. 36E-3 

0. 256 

6 

yi. is 

0. 35 

0. 02E-5 

0. 61E-4 

0, 18E-3 

0. 201 

— 

— 

0. 18E-3 

0.434 

i 

67. IS 
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indicate fairly close agreement with the accuracy estimates of 5 meters in 

range and 0. 006 degrees in azimuth and elevation. The cr value of 0. 0165 

VE 

degrees in Table V for Radar 7. 18 probably results from roughness in the 
data associated with the low elevation coverage (E° ^ 5. 33°) . 


CONCLUSIONS 


The TEMS Multiple Regression Analysis Method provides for a com- 
prehensive evaluation of systematic errors in measurements obtained from 
various tracking systems. The error model equations used in the develop- 
ment of the method are for C-band radar tracking systems. It should be 
noted, however, that the development for application to tracking systems 
other than radars is analogous to that presented herein. Truncated tracker 
error models for representing the systematic errors are established using 
the TEMS method in conjunction with a stepwise regression procedure. The 
stepwise procedure involves examining at every step the variables incorporated 
into the error model in previous steps. A specific variable is deleted from 
or entered into the model by using the Gaussian Elimination Method for solving 
the linear system of normal equations. The procedure shows considerable 
promise in solving the TEMS error model construction problem. 

The approach given by guidelines (1) - (3) in the last section for 
obtaining truncated error models to describe the systematic errors has 
generally resulted in acceptable models for the AS-501 first and second burn 
data. Guideline (2) uses the stepwise regression procedure to determine the 
variables for consideration in the truncated TEMS error models. It is noted 
that the average random errors remaining in the first and second burn residuals 
for the truncated error models indicate fairly close agreement with the input 
accuracy estimates of 5 meters in range, and 0. 006 degrees in azimuth and 
elevation. The stepwise results on the AS-501 data indicate several cases 
where the or curve fit value is not improved significantly by the introduction of 

additional variables into the regression. It appears that a rather critical 
examination of results from applications of the stepwise regression procedure 
is required in order to obtain meaningful and useful information for input to 
the TEMS program. 
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The scope of future work will be concerned with the reduction and 
analysis of tracking data on each Saturn IB and Saturn V flight tests. A 
comprehensive utilization of the stepwise regression analysis results will be 
stressed in this work. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, June 20, 1968 
933-50-04-00-62 
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APPENDIX A 

THE C-BAND RADAR TRACKING SYSTEM ERROR MODELS 
The Basic Error Models 


The basic radar error model equations [2] in range, azimuth, and 
elevation are given by: 

Range 

AR = C 0 + Ci R + C 2 R + C 3 t + C 4 (-. 022 cosec E) 


+ Cs (f) + c« (I) - c, (f ) 

( A-l) 

Azimuth 


AA — D 0 + Dj A + D 3 A + D5 tan E + D 3 sec E + tan E sin A 


• *r*v -m * . ^ t sin A cos A \ . ^ ( sin A cos A 

+ D 8 tan E cos A + D 9 ( — ) + D 10 ( — 


( A-2) 


+ D 41 A sec E 


Elevation 


AE - Fq + Fj E + F 3 E + F 5 (-sin A) + F e cos A 


+ F 7 


' . 022 
iR sin E 


- 10 6 ) cotan E 


+ Fq 


< -X tan E 
1 R 2 


( A-3) 


+ F 


10 


/ -Y tan E \ „ / 

y R 2 j + 


cos E 
R 


+ F 12 E cos E 
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The error model terms appearing in equations (A-l) , (A-2) , and 
(A-3) are subject to specific physical interpretations. These interpretations 
are listed below: 

In AR 

(1) C 0 - Zero set (bias) error which arises because of inaccuracies 
in calibration. 

(2) C t R - Range is obtained by multiplying the measured quantity by 
a scale factor. This is error due to error in the scale factor. 

(3) C 2 R - Error due to timing delays and is essentially the error in 
station time clock. 

( 4) C 3 t - Rate bias error due to linear drift of oscillator frequency 
with time. 

(5) C 4 ( -0.022 cosee E) - Error in correction made for tropospheric 
refraction. 

(6) C S (X/R), C 6 ( Y/R) , C 7 (Z/R) - Survey errors in the radar site 
location. 

In AA 

(1) D 0 - Same interpretation as for C 0 . 

• • 

(2) D 4 A - Same interpretation for C 2 R. 

(3) D 3 A* - Error due to dynamic lag in azimuth serve system. 

(4) D g tan E - Error which accounts for the nonpar allellism between 
the elevation shaft and the plane of the azimuth table. 

( 5) D e sec E - Error due to nonperpendicularity of the R-F axis 
and elevation shaft. 

( 6) D 7 tan E sin A, D 8 tan E cos A - Errors in azimuth due to 
azimuth plane tilt about orthogonal axes in the azimuth plane. 

(7) D 9 (sin A cos A/X) , D 10 (-sin A cos A/Y) - Errors in azimuth 
due to north and east components of survey error. 


64 



In AE 


( 8) Djj A sec E - Azimuth velocity lag. 

(1) F 0 - Same interpretation as for C 0 . 

(2) F 7 E - Same interpretation as for C 2 R. 

(3) F 3 E - Error due to dynamic lag in elevation servo system. 

(4) F 5 ( -sin A) , F e cos A - Error in elevation caused by azimuth 
plane tilt - same interpretation as for D 7 , D 8 . 

(5) F 7 ( ° 2 1 -10 -6 ) cotan E - Error in correction made for 

' 1_\ R sin E / J 

tropospheric refraction. 

(6) F 9 ( -Xtan E/R 2 ) , F 10 ( -Y tan E/R 2 ) , Fj^cosE/R), - Survey 
errors in radar site location. 

(7) F 12 E cos E - Elevation velocity lag. 

Derivation of Survey Error Terms in Range Error Model 


The error model equations contain terms which can, theoretically, be 
determined on the basis of physically established relationships. For example, 
consider the survey error coefficients C 5 , C 6 , and C 7 in the range error model 
equation. Let C 5 denote the error in the X direction, C 6 the error in Y direc- 
tion, and C 7 the error in the Z direction. Let the correctly measured range 
be given by: 

R = (X 2 + Y 2 + Z 2 )® . (A-4) 


If the radar site location is displaced by the distances AX, AY, and AZ in the 
respective coordinates, then the range from the site location (X + AX, Y + AY, 
Z + AZ) is given by: 

_i 

R = [ (X+ AX) 2 + (Y+ AY) 2 + (Z + AZ) 2 ] 2 . (A-5) 


65 


Expanding ( A - 5) : 

R = [X 2 + Y 2 + Z 2 + 2 (X AX+ Y AY + Z AZ) 

1 

+ (AX 2 + AY 2 + AZ 2 )] 2 . (A-6) 

For the application herein, the last three terms in equation (A-6) are small 
compared to the terms (X 2 + Y 2 + Z 2 ) and 2(X AX + Y AY + Z AZ) . Thus, 
these terms can be neglected and equation (A-6) written as: 

R = [X 2 + Y 2 + Z 2 + 2 (X AX + Y AY + Z AZ)]^ . (A-7) 

Expanding the right side of equation (A-7) in the binomial series and retaining 
the first two terms yields: 


B = ( X 2 + Y 2 + ZV * - XAX * YAY * Z , AZ . 

(X 2 + Y 2 + Z 2 ) 2 (A-8) 

The error in R is then given by: 

AR = R ‘ R m * (A " 9) 

Substituting equations (A-8) and (A-4) in(A-9): 



where 

AX = C 5 
AY = C 6 
AZ = C 7 


66 



APPENDIX B 


COORDINATE SYSTEM TRANSFORMATIONS 


The reference trajectory used as the standard in the TEMS method is 
in an earth-fixed plumbline coordinate system (X , Y , Z ) with origin at 

G G 6 

the launch site. It is required to transform the (X , Y , Z ) reference 

n e’ e e 

Y Y Y 

trajectory data into radar reference tracking parameters ( R , A , E ) . 

The radar measured tracking parameters can then be compared with the 
radar reference tracking parameters to establish the tracking errors for the 
particular system under consideration. There are two transformations required 
before this comparison can be made. These are: 

(1) Transformation of earth-fixed plumbline coordinates (X , Y , Z ) 

0 0 6 

of vehicle with origin at the launch site to earth-fixed plumbline coordinates 
(X , Y , Z ) of vehicle with origin at the tracking site. 

6S 6S 6S 

(2) Transformation of earth-fixed plumbline coordinate (X , Y , Z ) 

es es es 
r r r 

of vehicle with origin at tracking site to spherical coordinates ( R , A , E ) of 
vehicle with origin at tracking site. 

The first of these transformations for the i-th observation is given by 
the following equation: 


K T $ T X T A L # L K L^ X e + r OL ^ 

i 


where 



(B-i) 


(B-3) 
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I- 


OL 


-r sin. B cos K 


-r_ sin B m cos 

L Li L 


T T T 

r cos B t + h 

(B-4) r__ = 

r_ cos B_ + h,_ 

L L L 

OT 

T T T 

r sin B sin K 

_ L Li Ju 


r T sin B t sin 1^ 


(B-5) 



^sin 

k l 

0 

cos K 

Ij 



1 

O 

o 

1 

ttl 

ii 

0 


1 

0 

( B-6) 

II 

0 -cos -sin $ 

Li JL 


—cos 

k l 

0 

sin K 

Li , 



0 sin* -cos 4> 

■ Li Ju — 


( B-7) 



sin A 

Li 

-cos A 

Li 

o - 



sin K t 

0 

cos K t 

X L = 

cos A 

Li 

sin A 

L 

0 

( B-8) 

g T = 

0 

1 

0 


0 

0 

1 J 



-cos K t 

0 

sin K t _ 


(B-9) 


<3? = 

T 


0 -cos $ T -sin $ T 
0 sin $ T -cos $ TJ 


( B-10) A m = 


sin A t -cos A t 0 
cos A t sin A t 0 


(B-li) 


The equations required for calculating r r ^, B^, and B T appearing in equations 
{ B-4) and (B-5) are given by: 


B = $ - # 

L L L 


( B-12) 
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( B- 13 ) 


r = &/(cos 2 \S r T + — — sin 2 ^ ) 2 

j-, li g 2 i- 1 


B = $ - ^ 

T T T 


r m = ^/(cos 2 ^_ + — sin 2 ) 2 


£2 

*L = ta "‘‘ ~W ^ $ L 


* T = tan -1 -p- tan $ T . 


The spherical coordinate transformation for the i-th observation is given by: 


R. = (X 2 + Y 2 + Z 2 ) 

1 es. es. es. 

111 


A r = [tan -1 (Z /X )] + K , 0 ^ Af ^ 360 ° 
1 es. es. T’ 1 

1 1 


E. = tan " 1 
1 


(X 2 + Z 2 ) 2 

es. es. 

1 1 . 
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An additional transformation is required to determine the (X, Y, Z) 
position of the vehicle in an earth-fixed ephemeris coordinate system with 
origin at the tracking site. This transformation is given by: 


X. 

1 


X T $ T K T X es. 

1 


(B-21) 


where 



X. 

1 



( B-22) 


and X is obtained from the transformation given by equation ( B-l) . 
es. 

l 
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APPENDIX C 


A STEPWISE REGRESSION EXAMPLE ILLUSTRATING THE 
GAUSSIAN ELIMINATION METHOD 


As pointed out earlier, the Gaussian Elimination Method is ideally suited 
for a stepwise regression analysis since it actually enters variables into the re- 
gression equation one at a time. A significant variable entered into the 
regression equation at an early step may, after several other variables are 
entered, be found to be insignificant. This insignificant variable is then deleted 
from the regression equation before another variable is considered for entry. 
Thus, the final regression equation includes only the significant variables. 

The complete method is best illustrated by taking a specific example and per- 
forming the step-by-step operations. Assume there are p = 4 independent 
variables to be considered for entry into the regression equation. The aug- 
mented matrix similar to equation (117) is: 


r il 

r i2 

r i3 

r i4 

r iY 


C i2 

C i3 

C i4 

~11 

r i2 

r i3 

r i4 

r iY 

1 

0 

0 

0 


r 21 

r 22 

r 23 

r 24 

r 2Y 

0 

1 

0 

0 


00 

r 32 

r 33 

r 34 

r 3Y 

0 

0 

1 

0 

(C 

r 41 

r 42 

r 43 

r 44 

r 4Y 

0 

0 

0 

1 


r Yl 

r Y2 

r Y3 

r Y4 

r YY 

0 

0 

0 

0 



Step 1 

The first step consists of entering into the regression the variable most 
highly correlated with the response. Thus the variable corresponding to the 

largest of r^ , r^ , r^ , and r^ is selected for entry. Assume that r^ 

1Y Zjl oY 4Y 4Y 

is the largest. Then the nondiagonal elements r^, r^, r g4> and r of 
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the fourth column are zeroed out. This is accomplished by dividing the fourth 
row by r^; multiplying the new fourth row by ( -r ) and adding the result 

to the first row; multiplying the new fourth row by ( -r ) and adding the result 

to the second row; multiplying the new fourth row by (-r ) and adding the 

result to the third row; and finally, multiplying the new fourth row by ( -r ) 
and adding the result' to the fifth row : 


Jli 

r' 

X i2 

r* 

13 

V\ 

i4 

r’ 

iY 

c Ji 

c f 

i2 

c' 

i3 

c! 

i4 











r il r i4 r 41 
r 44 

r i2~ r i4 r 42 

r 44 

r i3 r i4f43 
r 44 

0 

r iY r i4 r 4Y 
r 44 

1 

0 

0 

-r /r 
14 44 

r 2l“ r 24 r 41 

r 44 

r 22" r 24 r 42 

r 44 

r 23~ r 24 r 43 
r 44 

0 

r 2Y _r 24 r 4Y 

r 44 

0 

1 

0 

-r /r 
24^44 

r 3i~ r 34 r 41 

r 44 

r 32 _r 34 r 42 

r 44 

r 33“ r 34 r 43 

r 44 

0 

r 3Y _r 34 r 4Y 

r 44 

0 

0 

1 

-r /r 
34 7 44 

r /r 
4f 44 

r /r 
42 44 

r /r 
43 44 

1 

P 4Y /r 44 

0 

0 

0 

l/r 

44 

r Yl"“ r Y4 r 41 

r 44 

r Y2 _r Y4 r 42 

r 44 

r Y3~ r Y4 r 43 

r 44 

0 

r YY~ r Y4 4 4Y 

r 44 

0 

0 

0 

- r Y/ r 44 


( C-2) 

At this step the fourth element in the fifth column is the solution a 4 in the 
normal equation: 


oi r = r 
4 44 4Y * 


( C-3) 


The fourth element in the fourth row of the matrix c' to the right of the 
line in equation (C-2) contains the inverse of r^. Of the three variables not 

in the regression the variable next considered for entry is the one whose 
partial correlation with the response is highest; that is, the variable corresponding 
2 2 2 

to the largest of p„ T ,, A , These would normally be determined 

as indicated earlier in this report. It can be verified, however, that ^ 

(i = 1,2,3) can be obtained using elements of equation (C-2) by the equation: 
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[ 


P iY.4 


(r f ) 
' iY ; 


11 YY 


( C-4) 


Where i = 1, 2, 3 and r' refers to the i-th diagonal element, r' refers to the 

ii iY 

element in the i-th row, fifth column and r' refers to the element in the fifth 

* * 2 

row, fifth column of equation (C-2) . Assume ^ is found to be the largest. 

Then the test for entering the variable is given by: 


PtY. 4 < n - d > 


1(IN) 


( C-5) 


1 -p 


1Y.4 


If this value is less than or equal to the appropriate table value, then 

the variable is not entered and the process is terminated. 


Step 2 


If, however, it is found that is greater than the table value, then 

the variable Z^ is entered. This is done analogous to the manner in which 
was entered in the first step. That is,, the nondiagonal elements r^, r^, r^, 
and r' of the first column in equation ( C-2) are zeroed out and the following 
new matrix is obtained: 


A" = 


r il 

r i2 

r i3 

r i4 

r iY 

cV 

U 

C i2 

c" 

13 

C i4 

1 

r 'l2 /r il 

r i 3 /r 'n 

0 

r iY /r il 

1/r ii 

0 

0 

c 1 / r 1 
14 11 

0 

r 22~ r 2i r i2 

r 23“ r 2i r l3 

0 

r 2Y _r 2i r iY 

- r 21 /r ll 

1 

0 

C 24" r 21 C 14 


r ii 

r ii 


r 'n 




r ii 

0 

r 32‘ r 3i r i2 

r 33" r 3i r i3 

0 

r* -r' r' 
3Y 31 1Y 

-r' /r' 
3i' 11 

0 

1 

C 34" r 31 C 14 


r ii 

r W 


r 1 
11 




r \i 

0 

r 42' r 4i r i2 

r 43' r 4i r i3 

1 

r 4Y" r 4i r iY 

-r’ /r 1 
41 11 

0 

0 

C 44 _r 41 C l4 


r il 

r W 


r \i 




r il 

0 

r» -r 1 r 1 

Y2 Y1 12 

r' -r 1 r T 
Y3 Y 1 13 

0 

r YY’ r Yi r iY 

-r 1 /r 1 
Y1 11 

0 

0 

C Y4“ r Yl C 44 


r' 

11 

r' 

11 


r il 




r il 


( C-6) 
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At this step the first and fourth elements in the fifth column are the 
solutions o'j and a ■ , respectively, in the normal equations: 



r ii 

r i4 


01 

1 


r 41 



OL 

4 

The elements in the il, 14, 
of the inverse matrix: 


" r il 

l-*» 

^ 1 

-1 



r 41 

r 44 


• 


•1Y 


4Y 


(C-7) 


Before the variables corresponding to a and a are considered for 

A u 

entry into the regression, the variables in the regression are tested for 
significance to determine if one should be deleted. This test for deletion is 
always made before a variable is added to the regression. The two variables 
in the regression corresponding to the coefficients a ^ and a ^ are considered 

for deletion by computing and F 4 ( G UT) from: 


F i( OUT) 


o>2 (n-d) 
^-<14 > 


(r^Hn-d) 

c*'. r" 
ii YY 


( C-8) 


where i = 1, 4. The smallest of F^ OUT ) and F 4 ^ OUT ^is tested for significance 

by comparison with the appropriate F table value. Assume the smallest is 
found to be greater than the appropriate table value. Then both Z^ and Z^ 
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remain in the regression and one of the two remaining variables is considered 
for entry into the regression. The variable considered for entry is the one 
whose partial correlation with the response is highest; that is, the variable 

2 2 

corresponding to the largest of p and p . It can be verified that 

2 ^ x • 14 oi. 14 

p. , (i= 2,3) can be obtained using elements of equation (C-6) by the 

i X • JL jl 

equation: 


2 

p iY. 14 


(r" ) 
iY 


-pTT -pH 

ii YY 


(C-9) 


where i = 2,3, and r" elements are the appropriate elements to the left of 

2 

the line in equation (C-6). Assume p is found to be the largest. Then 

£ X • It: 

the test for entering Z is given by: 


F 2(IN) 


P 2Y. 14 (n ~ d) 
. 2 
” P 2Y. 14 


( C-10) 


Step 3 


Assume F^ ^ is indicated to be significant from the comparison with 
the appropriate F table value. Then the variable Z corresponding to a is 
entered. The nondiagonal elements r" r" , r" , and r" of the second 

o£ xZ 

column in equation (C-6) are zeroed out and the following new matrix is 
obtained: 
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r'i» 

n 

J.I ft 

i2 

r tl» 

13 

r tii 

i4 

{.III 

iY 

c' Tf 

il 

pin 

C i2 

c ,,f 

13 


1 

0 

r i3" r i2 r 23 

0 

r lY- r i2 r 2Y 

_ r " c" 
11 12 21 

' r l2 /r 22 

0 

C 14~ r i2 C 24 



r 22 


r 22 

r" 

22 



r 22 

0 

1 

r *« / r »» 
23 22 

0 

r fi / r ir 

2Y 22 

c" / r" 
21 22 

1/C 22 

0 

C 24 /r 22 

0 

0 

r 33~ r 32 r 23 

0 

I, 3Y“ r 32 r 2Y 

C 3l" r 32 C 21 

" r 32 /r 22 

1 

C 34 _r 32 C 24 



r 22 


r 22 

r 22 



r 22 

0 

0 

r 43~ r 42 r 23 

1 

r 4Y _r 42 r 2Y 

C 4r r 42 C 21 

-r" /r" 
42 22 

0 

C 44‘ r 42 C 24 



r 22 


r 22 

r 22 



r 22 

0 

0 

r Y3“ r Y2 r 23 

0 

r ’YY" r Y2 r 2Y 

°Yr r Y2 C 21 

-r" / t" 
Y2 22 

0 

C Y4 _r Y2 C 24 



r 22 


r 22 

r 22 



r 22 

— 











(C-il) 


At this step the first, second, and fourth elements in the fifth column 
solutions a ay and ay respectively, in the normal equations: 


r , . r , r . 


a 


r . 

11 12 14 


1 


1Y 

r „ r r , 


QL 


r 

21 22 24 


2 


2Y 

v . r r . 


a 


r 

41 42 44 


4 


u 4Yj 


Elements of the inverse matrix 


r, . 

r , 

r , 

11 

12 

14 

r , 

r 

r 

21 

22 

24 

r „ 

r 

r 

41 

42 

44 


are the 


(C-12) 
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are contained in the appropriate positions of the c ,n matrix in equation (C-ll). 
The three variables in the regression are considered for deletion by computing 


F l(OUT) ’ 


F 2(OUT)’ and F 4(OUT) from: 


OL . 


(n-d) 


F i( OUT) 


°U (1 - B Y.124> 


<rj”) 2 (n-d) 


r» f t r y» f t f 

C ii r YY 


( C-13) 


where i - 1,2,4. The smallest (i=l,2,4) is tested for significance 

by comparison with the appropriate F table value. If the smallest is found to 
be greater than the appropriate F table value then all three variables correspond- 
ing to a a ^ would remain in the regression. The remaining variable 

would then be tested for entry by first computing: 


2 

P 3Y. 124 


( r ,M ) ‘ 
3Y 


33 YY 


( C-14) 


The r ,M elements are the appropriate elements in the matrix to the left of the 
line in equation (C-ll) . The test for entering Z is then given by: 

O 

„ _ 4r. 124 <n ' d) 

F 3( IN) 2 . (C-15) 

" P 3Y. 124 


If FgQ^s is found to be insignificant, 
final regression equation contains ct^ , 


then the process is terminated and the 

o' and ci . 

2 4 
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Step 4 


If, however, ^ is found to be significant, then the variable 
corresponding to a is entered by zeroing the nondiagonal elements r' n , 

o 13 

r"' r"' and r'" of the third column in equation (C-ll) . The following 

Zo 4o jlo 

new matrix is then obtained: 


IV 

IV 

r i2 

IV 

r i3 

IV 

r i4 

IV 

r iY 

~IV 

C il 

~IV 

IV 

13 

-IV 

_fl4_ 

1 

0 

0 

0 

-1 11 _ r «11 r 1 1 1 
1Y 13 3Y 

c n’ r i3 c 3i 

pi M_ r lltptH 

C 12 r i3 32 

13 x 33 

S 14' r l3^ 





r 33 

r"' 

33 

r 33 


r 33 

0 

i 

0 

0 

T»1 * * _T' M 1 * 

2V 23 3Y 

pit l_wt 1 1n HI 

21 23 32 

C 22" r 23 C 32 

' r 23' /r 33 

C nf_ r it i c nt 

24 23 34 





r '" 

33 

r .». 

33 

r in 

33 


r 1 H 

33 

0 

0 

1 

0 

r 1M/ r 1.1 

oY 33 

3/33 

C 32 /r 33 

l/r 33 

°34 /r 33 

0 

0 

0 

1 

r iit _ r m r m 

4Y 43 3Y 

p II 1 r» T » Ini 1 1 

C 41 r 43 C 31 

C 42" r 43 C 32 

-r’l’/r"! 
43 33 

c'i ,, - r , ”c ,M 
44 43 32 





r 33 

r lll 

33 

r 33 


33 

0 

0 

0 

0 

r m _ r m r i t» 

YY Y3 3Y 

C Yl _r Y3 C 31 

C Y2 _r Y3 C 32 

-m /r ,M 
Y3 7 33 

pill _ r 1 II pill 

Y4 Y3 32 





r m i 

33 

r tn 

33 

r 33 


r V 1 1 

33 


( C-16) 

At this step all the variables are in the regression and equation (C-16) is 
given by: 
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A IV = 


1 

0 

0 

0 

“l 

~IV 

C il 

~IV 

°12 

~IV 

°13 

~IV 

C 14 

0 

1 

0 

0 

a 2 

~IV 

C 21 

~IV 

C 22 

~IV 
C 23 

~IV 

°24 

0 

0 

1 

0 

Ol 

3 

~IV 
° 31 

~IV 

C 32 

~IV 

C 33 

~IV 

C 34 

0 

0 

0 

1 

“4 

~IV 

°41 

~IV 

°42 

~IV 

°43 

~rv 

°44 

0 

0 

0 

0 

( 1_R Y. 1234 ) 

~ a l 

~ a 2 

“ a 3 



(C-17) 


where: 


~IV 

~IV 

~IV 

~IV 

C 11 

°12 

C 13 

C 14 

~IV 

~IV 

~IV 

~IV 

C 2i 

C 22 

°23 

C 24 

~IV 

~IV 

~IV 

~IV 

C 31 

C 32 

C 33 

°34 

~IV 

~IV 

~IV 

~IV 

C 41 

°42 

°43 

°44 


r il 

r i2 

CO 

r i4 

r 21 

r 22 

r 23 

r 24 

r 31 

r 32 

r 33 

r 34 

r 41 

r 42 

r 43 

r 44 


( C-18) 


The four variables in the regression are considered for deletion by computing: 

2 


F. 


i(OUT) 


a (n - d) 


~IV , 2 

C ii 1 “ R Y. 1234 ) 


, IV ,2 . 
(r yy ) (n-d) 

~LV IV 
c. r 
ii YY 


( C-19) 
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where i = 1, 2, 3, 4. 


If the smallest F 


i(OUT) 


value is found to be greater 


than the appropriate F table value, then all four variables remain in the 
regression and the process is terminated. If, however, the smallest -F\( 

value is found to be insignificant, then the variable corresponding to this 


F. 


i(OUT) 


value is deleted from the regression by zeroing the appropriate 


elements in the column to the right of the line in equation (C-16) . The column 
selected corresponds to the particular variable being deleted. 


The deletion process will be illustrated by going back to equation (C-il). 
This equation contains the solution for the coefficients a a^, and a^. Assume 

^2(OUT) aS com P u * ;ec * by equation (C-13) is smaller than 

F (OUT) and is found to be insignificant. Then Z is deleted by zeroing the 
4 ~ * 

first, third, fourth, and fifth elements of the c"' column. This is accomplished 

’ ’ i2 

by dividing the entire second row by c"'; multiplying the new second row by 

( -c ,M ) and adding the result to the first row (zeros first element) ; multiplying 

the new second row by (-c'") and adding the result to the third row (zeros 

the third element) ; multiplying the new second row by ( -c^”) and adding the 

result to the fourth row (zeros the fourth element) ; and finally, multiplying the 
new second row by ( -c^) and adding the result to the fifth row ( zeros the 
fifth element) : 


V 

r. . 

u 

V 
1 i2 

V 

r i3 

V 

r i4 

V 

_^iY__ 

~v 

c.. 

li 

~V 

C i2 

-V 

°13 

~v 

C i4 

1 


*,» t tplll 

13 12 23 

0 

r iY“ C i2 r 2Y 

C il“ C 12 C 21 

0 

0 

C 14‘ C 12 C 24 



C 22 


C 22 

°22 



C> ' i 
22 

0 

1/? 22 

r 23 / "22 

0 

r IM /c ,M 
2 Y 22 

21 22 

1 

0 

c ,l, /c'" 
24 22 

0 

C 32 /C 22 

y t M_q n tr.1 M 

33 32 23 

0 

r m — c t i t r » » i 

3Y 32 2Y 

C 3l" C 32 C 21 

0 

1 

n f M f T f ft f 1 f 

34 32 24 



C 22 


C 22 

C 22 



C 22 

0 

C 42 X 22 

r 43' C 42 r 23 

1 

r 4Y _C 42 r 2Y 

C 4l’ C 42 C 2l 

0 

0 

C 44~ C 42 C 24 



°22 


TTt T t 
°22 

c'" 

22 



C 22 

0 

.pill /c' n 

°Y2 /C 22 

y t t t T 1 1 yt M 

Y3 Y2 23 

0 

r i 1 i _c n 'r'" 

YY Y2 2Y 

c ,M -C ,M C ,M 
Y1 Y2 21 

0 

0 

p t 1 f Iff p ? t f 

Y4 Y2 24 



pm 

C 22 


pi T » 
°22 

pm 

C 22 



yttf 

22 _ 


( C-20) 
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By substituting the elements c!” c\", c!" r!" and rV' of equation 

ll i2 i4 i3 lY 

(C-il) (where i = 1, 2 , 3, 4, Y) into equation (C-20) , the following equation 
is obtained: 



(C-21) 


This equation is seen to be equivalent to equation ( C-6) , which contains the 
solution for the variables and in the regression. 
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APPENDIX D 


TEMS IBM 7094 COMPUTER PROGRAM 
FOR C-BAND RADAR TRACKING SYSTEMS 


The various computational procedures in the TEMS method are 
implemented through the use of a highly flexible IBM 7094 computer program. 
Basically, the program consists of a main program and five subroutines. 

Control parameters are input by way of data cards in the main program to 
indicate the specific subroutines required for the reduction. Any combination 
of terms in the total error models can be retained in a given regression through 
the use of program control matrices. Trajectory data are processed independent 
of storage capacity, and bad data or discrete points can be deleted by the use of 
program control cards. A variable sampling rate can also be accepted. Addi- 
tional features of flexibility include variable plot scaling and double precision 
computations. 

One of the areas with considerable flexibility is in the selection of the 
specific error models for the reduction. Although constraints interrelate the 
equations for AR, AA, and AE, the results for a single equation regression 
can be obtained through the use of appropriate program control matrices. An 
18 x 18 control matrix AA(I, J) , consisting of l's and 0's, is used for selecting 
_ X ~ 

the elements of (B WB+W) to use in a specific computer run. An 18 x 1 control 

matrix BB(1, 1) is also used for selecting elements of (B^W N - B^WBC - We) 
to use. The use of these matrices enables any combination of terms in the 
total error models to be retained in a given regression. A summary of the 
basic logic and computations in the program is presented in the flow chart 
herein (Fig. D-l). Additional information on all aspects of the program 
development is given in Reference 1. 
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APPENDIX E 


UNIVAC 1 108 COMPUTER PROGRAM FOR APPLICATION OF 
THE STEPWISE REGRESSION ANALYSIS PROCEDURES 


This appendix contains a summary of the basic logic and computations 
in the UNIVAC 1108 MAIN and STEPRG programs (Fig. E-l) . The STEPRG 
program was developed by debugging and modifying a UNIVAC 1108 regression 
subroutine. It was made stepwise by incorporating into the subroutine the 
stepwise procedures presented eariler in this report and by appropriate use 
of the output. 
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INITIALIZATION FOR EACH STEP: 
0) b, =0, 1 = 1, 2 .... p 
(2) VOUT =-~ NOUT = 0 
YIN = 0 NIN = 0 

T=0 




FOR SMALLEST NEGATIVE V 


VOUT = V, 
NOUT =1 


N0 HAS I =p BEEN 

— — 1 DONE? 


°o = t - bl ^ '■■■ ' b P 

FOR p VARIABLES IN REGRESSION 


F UOUT1 =VOUT (df) 



DELETION | 

F OUT> I F i(OUT) | ^ 

K = NOUT 
df = df + 1 


r 0UT 1 H(OUT) 



CALCULATE 
NEW A(I,J) 
MATRIX 


DETERMINE THE COMPUTED 
RESPONSE: 


V V Y >™.^ 



I 


FIGURE E-i. (Concluded) 














APPENDIX F 


RESULTS FROM THE APOLLO-SATURN 501 
VEHICLE FLIGHT TEST 


This appendix presents a summary of the results from the Apollo- 
Saturn 501 Vehicle Flight Test launched on November 9, 1967. The Stepwise 
Regression Analysis results for the first and second burn data are presented 
in Tables F-I through F-VTI and Tables F-VHI through F-XI, respectively. 
Coefficient correlations for the truncated error models for the first and 
second burn data are given in Tables F-XII and F-Xin, respectively. 

In the figures which follow these tables (F-l through F-22) , the tracking 
errors for the various radars are represented by dots. The description of these 
errors as obtained from the TEMS least squares adjustment program is 
represented by the solid computed curves. 

The least squares residuals for the truncated error models presented 
in this appendix can be thought of as being composed of random errors and 
unmodeled systematic errors. A high random error content in the data may 
prevent a systematic error of comparable magnitude from being determined. 

The latter errors are those that can be attributed to uncertainties in the 
standard used in establishing the tracking errors, unknown systematic errors 
not absorbed by those that are modeled, or to geometry limitations. The 
presence of a significant unmodeled systematic error may prevent an adequate 
description of the data from being obtained. 
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TABLE F-I. RADAR 0. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 FIRST BURN DATA 


Variable 

No. 


r[C 2 | AR] 
r[ C 4 | AR] 
r[ C 5 ! AR] 
r[ C 6 | AR] 
r[ C 7 ! AR] 
r[ C 8 I AR] 


AR] = 0.76 
AR] = 0.94 
AR] = 0.16 
AR] = 0.96 
AR] =-0.33 
AR] =-0.96 
AR] = 0.87 


Step 

No. 

Variables in Regression 

a Y 

F Level 


Final 


Co. Cg, C 7 
Co, C 5 , C 7 , Cg 
Cq, Cg, C 7 , Cg, c 4 


Cq, Cg, C 7 , Cg 


c 4 , c 6 


Co* Cg, Cg, C 4 , Cg 


4459. 1 
301.3 
204. 1 
9.3 
4. 9 


r[C 2 

r[D 3 

r[D 5 

r[D e 

r[D7 

r[D 8 

r[C 5 

r[C e 


I AA] 

! AA] 

I AA] 

I AA] 

I AA] 
i AA] 

! AA] 

I AA] 


0.60 
0. 55 
0. 09 
0. 12 
-0.43 
0. 50 
-0.03 
=-0. 13 


D 0 , C 2 , Dg 
Do, C 2 , D 3 , D 7 
Do, C 2 , Dg, D 7 , Dg 
Dq, C 2 , D 7 , Dg 


Final . D 0 , C 2 , D 7 , D 8 , Cg 


0. 0089 
0. 0062 
0. 0061 
0. 0058 
0. 0058 


0. 0058 


191.5 

353.9 

13.6 

30.9 

- 0.2 


= - 0.22 
= 0. 15 
= -0. 31 
= -0.44 
= 0.01 
= 0.49 
= 0.06 
= -0. 35 


F 0 , C 5 0.0074 

F 0 , Cg, Cg 0. 0067 

Fo» Cg, Cg, C 2 0.0061 

Fo, Cg, Cg, C 2 , D 7 0.0060 

Fq, CE, Cg, C 2 , D 7 , Dg 0.0060 


Fq, Cg, C 2 , D 7 , Dg 0,0060 


104.5 
68.5 
74.9 
9. ‘ 
3. ( 










TABLE F-II. RADAR 19. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 FIRST BURN DATA 


CD 




TABLE F-IH. RADAR 7, 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 FIRST BURN DATA 


Variable 

No. 

Linear Correlation 

Step 

No. 

Variables in Regression 

a Y 

F Level 

1 

r[ Ci | AR] = -0. 36 

1 

Co, C 7 

3.35 

3089.3 

2 

r[ C 2 | AR] = 0. 15 

2 

C 0 , C 7 , C 6 

2.86 

110.7 

3 

r[ C 4 | AR] = 0. 77 ; 

3 

Co, C 7 , C 6> C 7 

1 . 16 

1493. 9 

4 

r[ C 5 | AR] = 0. 31 

4 

Co, C 7 , Cg, Cj, Cg 

1.15 

4.9 

5 

r[ C e 1 AR] = 0. 32 





6 

r[ C 7 1 AR] = 0. 96 





7 

r[ C 8 1 AR] = 0. 36 

Final 

Co, C 7 , Cg, Cj, Cg, C 2 

1.13 

15.0 

1 

r[C 2 1 AA] = 0.28 

1 

Do, D 7 

0. 0037 

76.6 

2 

r[ D 3 | AA] = - 0 . 44 





3 

r[ D 5 1 AA] = - 0 . 09 





4 

r[ D e | AA] = - 0 . 09 





5 

r[D 7 | AA] = 0.45 





6 

r[ D 8 1 AA] = -0. 06 




! 

7 

r[ C 5 1 AA] = -0. 05 




i 

8 

i r[ C 6 1 AA] = -0. 37 

Final 

Only 1 step. 



1 1 

' r[C 2 1 AE] = 0. 03 

1 1 

Fo> D 7 

0.0158 

22.9 

2 

r[ F 3 1 AE] = 0. 22 




i 

3 

j r[ D 8 1 AE] = 0. 01 




1 

4 

r[ D 7 1 AE] = -0. 27 





5 

i r[ C 4 1 AE] = 0. 07 




1 

j 

6 

r[C 5 1 AE] = 0.004 





7 

r[ C 6 1 AE] = 0. 24 





8 

r[ C 7 t AE] = -0. 16 







Final 

Only 1 step. 





TABLE F-IV. RADAR 3. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 FIRST BURN DATA 


r 

Variable 

No. 

1 

1 

Linear Correlation 

Step 

No. 

| 

Variables in Regression 

a Y 

F Level 

1 

r[ C 4 | AR] = 0.47 

1 

C 0 , C 5 

5.61 

1669. 0 

2 

r[C 2 | AR] = 0.73 

2 

C 0 , C 5 , Cj 

2.36 

1834.7 

3 

r[ C 4 | AR] = 0. 25 

3 

Co, c 5 , C 4 , c 7 

2. 03 

136. 1 

4 

r[ C 5 | AR] = 0. 90 

4 

Co, Cg, C 4 , C 7 , Cg 

1.93 

45.3 

5 

r[C 6 | AR] = 0.62 

5 

Co, Cg, C u C i, Cg, c 2 

1.87 

24.2 

6 

r[ C 7 | AR] = -0. 26 





7 

r[ C 8 | AR] = 0. 74 

Final 

Co, Cg, C 4 , C 7 , Cg, C 2 , Cg 

1.82 

22.8 

1 

r[ C 2 | AA] = 0. 62 

1 

Do, c 2 

0.0028 

243.7 

2 

r[ Dg | AA] = 0.43 

2 

Do, c 2 , Dg 

0. 0025 

117. 1 

3 

r[ D 5 | AA] = 0. 30 

3 

Do, C 2 , Dg, Dg 

0. 0024 

7. 9 

4 

r[ D e | AA] = 0. 32 





5 

r[ D 7 | AA] = -0. 60 





6 

r[D 8 | AA] = 0.41 





7 

r[ C 5 | AA] = 0. 07 





8 

r[ C 6 | AA] = 0. 51 

Final 

Do, C 2 , Dg, Dg, D 7 

j 0.0024 

35. 0 

1 

r[C 2 | AE] = “0. 11 

1 

F 0 , C 6 

0. 0053 

104. 1 

2 

r[ F 3 | AE] = 0. 44 

2 

Fo, Cg, Cg 

0.0051 

25.2 

3 

r[ Dg | AE] =-0. 08 

3 

Fo, Cg, Cg, D 7 

0. 0044 

133.9 

4 

r[ D 7 | AE] =-0. 42 

, 4 

F 0 , Cg, Cg, D 7 , C 4 



5 

r[C 4 | AE] = -0. 13 





6 

r[C 5 | AE] = -0. 01 





7 

r[C 6 | AE] = 0.46 





8 

r[C 7 [ AE] =-0. 23 







Final 

Fq, Cg, Cg, D 7 , C 4 , Fg 

0. 0043 

16.8 



CD 

4 ^ 


TABLE F^V. RADAR 67. 16 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 FIRST BURN DATA 


Variable 

No. 

Linear Correlation 

Step 

No. 

Variables in Regression 

°Y 

F Level 

1 

r[ Cj | AR] =-0.88 

1 

C 0 , C 7 

14. 66 

990.2 

2 

r[C 2 | AR] =-0.49 

2 

C 0 , C 1? C 6 

7.96 

687.8 

3 

r[ C 4 | AR] = 0. 83 

3 

Co> Ci, Cg, c 7 

3.91 

898.3 

4 

r[ C g | AR] =-0. 25 

4 

Cq, Cg, c 7 

3.91 

-1.6 

5 

r[C 6 | AR] =-0.78 

5 

Co, C 6 , C 7 , Cg 

2.76 

291.5 

6 

r[C 7 | AR] = 0.23 





7 

r[C 8 | AR] =-0.31 

Final 

Co, c 6 , C 7 , C 8 , c t 

2.73 

7.5 

1 

r[ C 2 | AA] = 0. 07 

1 

Do, D 3 

0. 0099 

1178.0 

2 

s> 

I 

II 

o 

CD 

O 

2 

Do, Dg, D 7 

0. 0078 

182. 9 

3 

r[ D 5 | AA] =-0. 10 

3 

Do, Dg, D 7 , D 5 

0. 0072 

47.7 

4 

r[ D e | AA] =-0. 10 





5 

r[ D 7 | AA] = 0. 84 





6 

r[ Dg | AA] = 0.25 





7 

r[ C g | AA] =-0.11 





8 

r[ C 6 | AA] = 0.03 

Final 

Do, D3, D 7 , Dg, Cg 

0. 0072 

3.8 

1 

r[C 2 | AE] = 0.46 

1 

F 0 . Fg . 

0. 0092 

299. 1 

2 

r[F 3 | AE] = 0.71 

2 

F 0 , Fg, D 8 • 

0.0050 

696.4 

3 

r[D 8 | AE] = 0.63 

3 

Fo, Fg, Dg, C 7 

0. 0045 

65. 0 

4 

r[ D 7 | AE] = 0. 05 





5 

r[ C 4 | AE] = 0. 02 





6 

r[ C 5 | AE] = 0. 70 





7 

r[C„ | AE] =-0.42 





8 

r[C 7 | AE] = 0.41 







Final 

Fq, Fg, Dg, C 7 , D 7 

0. 0044 

10.1 


/ 









TABLE F-VI. RADAR 67. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 FIRST BURN DATA 


Variable 

No. 

Linear Correlation 

Step 

No. 

1 

Variables in Regression 

°Y 

F Level 

i 1 

rfCj | AR] =-0.68 

1 

Co, C 6 

18.03 

1 394.5 

2 

r[C 2 | AR] =-0.32 

2 

Co, Cg, c 4 

3.61 

17061. 2 

3 

r[ C 4 | AR] = 0. 60 

3 

Cq, Cg, Cj, C 7 

2. 18 

511.4 

4 

r[C 5 | AR] =-0.004 

4 

Co, Cg, c 4 , c 7 , c 4 

2.04 

41.9 

1 5 

r[C 6 | AR] =-0.76 

5 

Cg, C 6 , c lf c 7 , C 4 , c 2 

1.95 

30.0 

1 6 

r[C 7 | AR] = 0.33 





7 

r[ C 8 | AR] = 0.26 

Final 

Co, C ls c 7 , c 4 , c 2 

1.95 

- 0.8 

1 

r[C 2 | AA] =-0. 14 

1 


0. 0043 

131.3 

2 

r[D 3 | AA] = 0.56 

2 

Do, D 3 , Dg 

0. 0042 

14.6 

3 

r[ D 5 |AA] = 0.15 





4 

r[D e ! AA] = 0.13 





■ 5 

r[ D 7 | AA] = 0. 38 





: 6 

r[ D g | AA] =-0.06 





7 

r[C 5 | AA] =-0.08 





8 

r[C 6 | AA] =-0.05 

Final 

Do, D 3 , Dg , D 7 

0. 0042 

5.6 

’ 1 

r[ C 2 | AE] = 0. 53 

1 

F 0 , C 2 

0. 0046 

118.2 

2 

r[ F 3 |AE] = 0.05 

1 

2 

Eo, C 2 , D 7 

0. 0044 

34. 0 

3 

r[ D g |AE] = 0.49 

3 

! F 0 , C 2 , D 7 , F 3 

0. 0038 

94.6 

4 

r[ D 7 |AE] =-0.47 

4 

Eo, C 2 , D 7 , F 3 , C 5 

0. 0035 

50.7 

5 

r[C 4 | AE] = 0.15 

5 

Eo, C 2 , D 7 , F 3 , Cg, C 7 

0. 0035 

12.4 

6 

r[C 5 | AE] = 0.19 





7 

r[C 6 | AE] = 0.26 





8 

r[ C 7 | AE] = 0. 21 







Final 

Eq, c 2 , d 7 , F 3 , Cg, C 7 , C 4 

0. 0034 

7.9 



TABLE F-VH. RADAR 1. 16 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 FIRST BURN DATA 


Variables in Regression 


Variable 

No. 


Linear Correlation 


r[Cj | 

AR] 

r[C 2 | 

AR] 

r[C 4 | 

AR] 

r[C 5 | 

AR] 

r[C 6 | 

AR] 

r[C 7 | 

AR] 

r[C 8 

1 AR] 

r[C 2 

1 AA] 

r[Dg 

1 AA] 

r[D 5 

1 AA] 

r[D e 

1 AA] 

r[D 7 

1 AA] 

r[D 8 

1 AA] 

r[C 5 

1 AA] 

r[C 6 

1 AA] 

r[C 2 

1 AE] 

r[F 3 

1 AE] 

r[D g 

1 AE] 

r[ D 7 

1 AE] 


r[ C 5 I AE] 


r[ C 7 I AE] 


C n. Cc. c 


^0* C 2 > U 7 * Dg 


Do, Cg, D 7 , Dg, C 
Fo 


Only 1 Step. 


0. 0126 
0.0106 
0. 0103 


0.0102 

0.0101 


F Level 


1190. 1 















TABLE F-Vffl. RADAR 91. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 SECOND BURN DATA 


Variable 

No. 

Linear Correlation 

— 

Step 

No. 

Variables in Regression 

ct y 

F Level 

1 

r[Cj |AR] =-0.97 

1 


4.12 

63402. 7 

2 

r[C 2 | AR] =-0.85 

2 

Co, C 8 , C 7 

2.61 

1024. 5 

3 

r[C 4 | AR] = 0.21 

3 

^2 

2.04 

438.3 

4 

r[C 5 | AR] =-0.88 

4 

^Oj ^4 

1.90 

99.4 

5 

r[ C 6 | AR] =-0.97 





6 

r[ C 7 | AR] = 0. 99 





7 

r[C 8 | AR] =-0.99 

Final 

C 8 , C 7 , C 2 , C 4 , Cj 

1.77 

105.8 

1 

r[C 2 | AA] = 0.14 

1 

Do, c 2 

0. 0037 

13.5 

2 

r[ D s | AA] = 0. 08 

2 

Do, C 2 , Dg 

0. 0037 

3. 9 

3 

r[D g | AA] =-0.01 

3 

Do, c 2 , Dg, D 6 

0. 0037 

4.7 

4 

r[ D g | AA] =0. 005 

4 

Do, c 2 , Dg, D e , D 7 

0. 0036 

7.8 

5 

r[D 7 | AA] =-0.14 

5 

Do, C 2 , Dg, D 7 

0. 0036 

-0.05 

6 

r[ D 8 | AA] = 0.11 





7 

r[C 5 | AA] = 0.09 





8 

r[ C 6 | AA] = 0. 002 

Final 

Do, C 2 , Dg, D 7 , D 3 

0.0035 

36.7 

1 

r[C 2 | AE] = 0.38 

1 

F 0» C 2 

0. 0055 

118.4 

2 

r[F 3 | AE] =-0.25 





3 

r[ D 8 | AE] = 0.37 





4 

r[ D 7 |AE] = 0.35 





5 

r[C 4 AE] =_0. 02 





6 

r[C 5 | AE] = 0.18 





7 

r[ C 6 AE] = 0.14 





8 

r[C 7 | AE] = 0.34 







Final 

Eq, C 2 , Fg 

0. 0054 

11.8 


TABLE F-IX. RADAR 3. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 SECOND BURN DATA 


Variable 


Linear Correlation 


Variables in Regression 


F Level 


AR] =-0.03 
AR] =-0.98 
AR] =-0.47 
AR] =-0.97 
AR] =-0.87 
AR] = 0.44 
AR] =-0.96 


c 6 , c 5 

^6, C 5 , C t 


6504.7 
160. 9 
214. 1 
221. 1 


AA] =-0.16 
AA] = 0.34 
AA] =-0.21 
AA] =-0.22 
AA] =-0.24 
AA] = -0. 14 
AA] =-0.01 
AA] = 0. 19 


D 0 , D 3 
Do, D 3 , D 7 
Do, D 3 , D 7 , D 5 
Dq, D 3 , D 7j D 5 , D e 


Dq, D 3 , D g , D e 


0. 0038 
0. 0037 
0. 0036 
0. 0035 


0. 0035 


r[C 2 
r[F 3 
r[D 8 
r[ D 7 
r[C 4 
r[C 5 
r[C 6 
r[C 7 


AE] =-0.08 
AE] = 0.53 
AE] = 0.08 
AE] =-0.58 
AE] =-0.51 
AE] =-0. 0C 
AE] = 0.52 
AE] =-0.67 


F 0 , C 7 


0. 0048 


267. 6 


Only 1 step 












TABLE F-X. RADAR 19. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 SECOND BURN DATA 


Variable 


S 

No. 

Linear Correlation 



Variables in Regression 


F Level 


rlCj 

[AR] =-0.63 

1 

C 0 , C 2 

8.29 

11057.6 

r[C 2 

| AR] =-0.98 

2 

Co, C 2 , c 8 

4.73 

1019. 6 

r[C 4 

| AR] = 0.27 





r[C 5 

| AR] =-0.97 





r[C 6 

| AR] = -0. 97 





r[C 7 

| AR] = 0.66 





r[C 8 

| AR] = -0. 98 

Final 

Co, C 2 , C 8 , C 7 

4.12 

153.8 

r[C 2 

| AA] =-0.48 

1 

^o, D 5 

0.0051 

185. 9 

r[D 3 

| AA] = 0.48 

2 

Do, Dg, D 3 

0. 0043 

198.9 

r[D 5 

| AA] = -0. 52 

3 

Do, Dg, D 3 , D 7 

0. 0042 

26.5 

r[D e 

| AA] =-0.50 

4 

Do, Dg, Dg, D 7 , D s 

0. 0042 

3.9 

r[D 7 

| AA] = -0. 37 

5 

Do, Dg, Dg, D 7 , D 8 , D e 

0.0039 

84.3 

r[D 8 

| AA] =-0.42 





r[C 5 

| AA] =-0.41 





r[C 6 

| AA] = -0. 00 

Final 

Do, D s , D 7 , D g , Dg 

0. 0039 

-1.5 

r[C 2 

| AE] =-0.18 

1 

F o, C 7 

0. 0068 

739.6 

r[F 3 

| AE] = 0.40 

2 

Fo, C 7 , C 4 

0. 0065 

44.3 

r[D 8 

| AE] =-0.35 

3 

Fo, C 7 , C 4 , D 7 

0. 0064 

26. 9 

r[D 7 

| AE] = -0. 75 

4 

Fo, C 4 , D 7 

0. 0064 

-1.4 

r[C 4 

| AE] =-0.76 

5 

Fo, C 4 , D 7 , F 3 

0. 0063 

14.7 

r[C 5 

| AE] = 0. 02 

6 

Fo, C 4 , D 7 , F 3 , D 8 

0. 0062 

13.4 

r[C 6 

| AE] = 0. 33 





r[C 7 

| AE] =-0.76 

Final 

Fflj C 4 , Fq, Do 

0. 0062 

-0.2 



u ? f± ~ o~ 0 









TABLE F-XI. RADAR 67. 18 STEPWISE REGRESSION ANALYSIS RESULTS FOR AS-501 SECOND BURN DATA 


Variable 


Step 


°Y 


No. 

Linear Correlation 

No. 

Variables in Regression 

F Level 



r[C t 

| AR] =-< 

r[Q 

| AR] =-< 

r[C 4 

I AR] = 

r[C 5 

| AR] =- 

r[C 6 

1 AR] =- 

r[C 7 

1 AR] = 

r[C 8 

| AR] =- 

r[C 2 

| AA] =- 

r[D 3 

1 AA] =- 

r[D 5 

1 AA] = 

r[D e 

1 AA] = 

r[D 7 

| AA] =- 

r[D 8 

| AA] =- 

r[C 5 

1 AA] = 

r t C 6 

I AA] =- 

r[C 2 

| AE] = 

r[F 3 

t AE] =- 

r[D 8 

| AE] = 

r[ D 7 

1 AE] = 

r[C 4 

1 AE] =- 

r[C 5 

| AE] =- 

r[C e 

| AE] = 

r[ C 7 

1 AE] - 


C 0 , C t 

8.73 

Co, C l5 C e 

4.89 

Co, C lt Cg, C 4 

4.76 

Co, C l5 Cg, c 4 , c 2 

4.58 

Co, c lt Cg, c 4 , c 2 , c 5 

4.23 

Co, Cg, c 4 , c 2 , c 5 

4.23 

Co, Cg, C 4 , C 2 , Cg, c 8 

4. 13 

Co, Cg, C 4 , C 2 , Cg, c 8 , c 7 

3.91 

All 

3. 85 

Do, C 2 

0. 0046 

Do, C 2 , Dg 

0. 0046 

Do, C 2 , Dg, D 7 

0. 0042 

Do, C 2 , Dg, D 7 , Dg 

0. 0041 

Do, C 2 , Dg, D 7 , Dg, Dg 

0. 0040 

Do, C 2 , Dg, D 7 , Dg, Dg, D 3 

0. 0040 

Fo, Dg 

0. 0057 

Fo, Dg, C 2 

0. 0054 

Fn. Do. Co. C, 

0. 0053 


28286. 9 
1889.6 











TABLE F-XEI. COEFFICIENT CORRELATIONS FOR THE 
TRUNCATED AS-501 FIRST BURN RADAR ERROR MODELS 


Co 

Cl 

C 2 

D 0 

d 3 

Dg 

F 0 

f 3 

I 1. 00 

-0.92 

0. 12 

0. 01 

0. 

-0. 02 

0. 

0. 

Cl 

1. 00 

-0. 16 

-0. 01 

0. 

0. 03 

0. 

0. 


c 2 

1. 00 

0. 05 

0. 01 

-0. 18 

0.01 

0. 



Do 

1. 00 

-0. 03 

0.25 

-0. 01 

-0. 01 




d 3 

1. 00 

-0. 12 

0. 01 

0. 01 





d 8 

1. 00 

-0. 04 

-0. 05 






F 0 

1. 00 

0. 03 







f 3 

1. 00 

C 0 

Ci 

Radar 67. 16 
C 2 Do 

d 3 

Dg 

F 0 

f 3 

1. 00 

-0. 88 

-0. 18 

-0. 01 

0. 

0. 05 

-0. 02 

0. 

Cl 

1. 00 

0.41 

0.02 

0. 01 

-0. 10 

0.04 

0. 01 


C2 

1. 00 

0. 05 

0. 02 

-0. 25 

0. 09 

0. 01 



D 0 

1. 00 

-0. 03 

0. 26 

-0. 10 

-0. 02 




d 3 

1. 00 

-0. 13 

0. 05 

0. 01 





Dg 

1. 00 

-0.39 

-0. 07 






F 0 

1. 00 

0. 06 








1. 00 


Radar 67. 18 
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TABLE F-XEI. (Concluded) 


Co 

c 2 

C 4 

Do 

d 3 

F 0 

f 3 

o 

o 

T-i 

- 0.37 

0.56 

0.04 

0 . 

0.02 

0 . 

c 2 

1 . 00 

0.41 

- 0. 12 

0 . 01 

0 . 02 

0 . 01 


C 4 

1 . 00 

- 0.05 

0 . 

0.04 

0 . 01 



D 0 

1.00 

0 . 04 

0 . 

0 . 




d 3 

1.00 

0 . 

0 . 





F 0 

1.00 

0.30 






*3 

1 . 00 


Radar 3. 18 


Co 

Cl 

c 2 

D 0 

d 3 

d 7 

F 0 

1 . 00 

0 . 19 

- 0.61 

- 0 . 11 

0.21 

0.27 

- 0.23 

Cl 

M- 

O 

O 

0 . 84 

- 0 . 15 

P 

to 

l " 

0.37 

- 0.31 


c 2 

1 . 00 

0. 18 

- 0.35 

- 0.44 

0 . 37 



D 0 

1 . 00 

- 0.56 

- 0 . 80 

0 . 73 




d 3 

O 

O 

• 

-T-i 

0.68 

- 0 . 62 





d 7 

O 

o 

- 0 . 91 






F 0 

1 . 00 


Radar 1.16 
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TABLE F-XIH. COEFFICIENT RADAR CORRELATIONS FOR THE 
TRUNCATED AS-501 SECOND BURN ERROR MODELS 



m 






















































TABLE F-XLH. (Concluded) 


C 0 

Cl 

C 2 

D 0 

d 3 

F 0 

f 3 

1. 00 

-0.67 

0. 01 

0. 

0. 

0. 

0. 

Cl 

O 

o 

1 

-0. 69 

0.04 

0. 

0. 

0. 


C 2 

1. 00 

-0.05 

0. 

0. 

0. 01 



D 0 

1.00 

0.35 

0. 

0. 




D 3 

1. 00 

0. 

0. 





F 0 

1.00 

0.44 






f 3 

1. 00 


Radar 91. 18 
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FIGURE F-4. RADAR 19.18 RANGE, AZIMUTH, AND ELEVATION 
ERRORS ON AS-501 FIRST BURN DATA 












TIME , SEC 


FIGURE P-7. RADAR 7. 18 RESIDUALS ON AS-501 
FIRST BURN DATA 
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FIGURE F-8. RADAR 7. 18 RANGE, AZIMUTH, AND ELEVATION 
ERRORS ON AS-501 FIRST BURN DATA 
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FIGURE F-10. RADAR 67. 16 RANGE, AZIMUTH, AND ELEVATION 
ERRORS ON AS-501 FIRST BURN DATA 
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FIGURE F-12. RADAR 67. 18 RANGE, AZIMUTH, AND ELEVATION 
ERRORS ON AS-501 FIRST BURN DATA 













FIGURE F-16. RADAR 19.18 RANGE, AZIMUTH, AND ELEVATION 
ERRORS ON AS-501 SECOND BURN DATA 





FIGURE F-17. RADAR 3. 18 RESIDUALS ON AS-501 
SECOND BURN DATA 






FIGURE F-18. RADAR 3. 18 RANGE, AZIMUTH, AND ELEVATION 
ERRORS ON AS-501 SECOND BURN DATA 
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FIGURE F-19. RADAR 19. 18 RESIDUALS ON AS-501 
SECOND BURN DATA 





FIGURE F-20 RADAR 1 9 . i 8 RANGE, AZIMUTH, AND ELEVATIf 
ERRORS ON AS-501 SECOND BURN DATA C 







FIGURE F-22. RADAR 67. 18 RANGE, AZIMUTH, AND ELEVATION 
ERRORS ON AS-501 SECOND BURN DATA 
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